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ABSTRACT 


In  Part  I of  this  report,  the  performance  of  a CENSORING  system 
is  examined  from  the  point  of  view  of  determining  the  probability,  P;#nsor 
that  at  any  instant  of  time  the  random  wavefront  distortion  over  a circular 
aperture  of  diameter  D will  be  small  enough  to  allow  a nearly  diffraction- 
limited  image  to  be  formed.  It  is  pointed  out  that  the  effective  wavefront 
distortion  for  this  calculation  is  not  the  total  distortion  function,  Ofr)  > 
but  rather  :c(r:L>)  which  represents  the  distortion  after  subtraction  of 
the  instantaneous  average  phase  and  tilt  over  the  aperture.  The  problem 
of  calculating  the  probability  is  related  to  the  probabilities  for  the  value 
of  the  various  random  coefficients  3b  of  the  decomposition  of  tp(r;D)  into 
a series  pr  fB(r;D)  , where  the  functions  f„  are  orthonormal  functions 
over  the  aperture  chosen  so  as  to  make  the  various  Sn  statistically  inde- 
pendent. (It  is  noted  that  the  3„  are  gaussianly  distributed,  since  3 is 
a gaussian  random  function.  ) The  appropriate  Karhunen- Loeve  homogeneous 
integral  equation  is  developed  to  allow  the  fB  to  be  obtained  as  eigenfunctions. 
The  variance  of  the  3,  are  seen  to  be  the  corresponding  eigenvalues.  It  is 
then  shown  how  the  probability  PCeBtor  can  calculated  as  a multi-dimen- 
sional integral  over  the  product  of  gaussian  distribution  with  the  variances 
corresponding  to  those  of  the  3n  . The  numerical  solution  of  the  Karhunen- 
Loeve  homogeneous  integral  equation  to  obtain  the  eigenvalues  (and  the  eigen- 
functions) and  the  numerical  evaluation  of  the  multi -dimensional  integral 
which  finally  gives  Pc„r.«or  are  taken  UP  'n  Part  II  this  report. 

In  Part  II  of  this  report,  the  eigenvalues  and  the  eigenfunctions  for 
average  tilt  and  average  phase  suppressed  wavefront  distortion  on  a circular 
aperture  are  developed.  (The  eigenvalues  and  eigenfunctions  without  tilt 
distortion  suppressed  are  also  developed.  ) Using  the  eigenvalue  set,  the 
CENSORING  system  probability,  PctNSOI1  , of  getting  a short  exposure 
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image  with  less  than  one  radian  squared  distortion  averaged  over  the 
aperture  is  evaluated  as  a multi-dimensional  integral,  evaluated  by  Monte 
Carlo  techniques.  The  results  are  found  lo  have  the  form 

Pcensor  * 5.6  exp  0.1557  (D/rQ)2  ] 

The  exponential  dependence  on  aperture  area  is  in  agreement  with  an 
earlier  conjecture  by  Hufnagel  (though  with  significantly  different  coef- 
ficients than  were  suggested  by  Hufnagel). 

It  is  noted  that  this  exponential  dependence  on  aperture  diameter 
makes  the  proper  selection  of  D/r0  a very  critical  aspect  of  a CENSOR- 
ING experiment.  It  is  pointed  out  that  if  independent  samples  of  wave- 
front  distortion  are  obtained  every  10  msec,  and  D/r0  = 15  (corres- 
ponding to  D = 1.  5 m , r0  = 0.  1 m ),  then  we  would  have  to  wait  about 
800  million  hours  for  a good  picture.  If  D/r0  is  reduced  to  a value  of 
10  (corresponding  to  D = 1.  0 m , r0  = 0.  1 m ),  the  waiting  time  is 
reduced  to  about  2.  5 hours,  while  if  D/r0  is  reduced  to  7 (corres- 
ponding to  D = 0.  7 m , r0  = 0.  1 m ),  the  waiting  time  drops  to  only 
3.  5 seconds. 

The  report  has  been  bound  in  two  volumes.  The  main  volume 
presents  the  derivation  and  principal  results.  The  addendum  volume 
presents  the  more  voluminous  tables  of  intermediate  results,  particularly 
those  that  may  be  of  use  in  working  other  problems. 


PART  I 


Forma!  Theory 
o f 

CENSORING  Systems  Opera 


t i o n 


Introduc!  ,on 


The  concept  of  CENSORING  as  a method  of  obtaining  high  resolution 
images  through  atmospheric  turbulence  is  based  on  the  assumption  that  of 
all  possible  forms  that  random  wavefront  distortion  will  assume  during  some 
period  of  time,  there  is  a finite  probability  that  at  some  instant  the  random 
pattern  wdl  very  nearly  resemble  a plane  wave.  At  that  instant,  a nearly 
diffraction-limited  image  can  be  obtained.  A CENSORING  system  would  be 
able  to  recognize  this  condition  quickly  enough  to  allow  a photograph  to  be 
taken  just  then,  while  preventing  exposures  at  times  of  more  normal  dis- 
tortion. 

The  simplicity  of  the  CENSORING  concept  makes  it  seem  quite 
attractive,  vut  on  the  other  hand  it  is  entirely  dependent  on  random  occur- 
rences for  its  operation.  It  is  therefore  critical  that  we  understand  the 
probabilities  involved  and  be  able  to  estimate  the  time  we  will  have  to  wait, 
on  ar  average,  before  the  CENSORING  system  will  provide  us  with  a picture. 
It  appears  likely  that  the  probability  that  at  any  instant  of  time  the  distorted 
wavefront  will  be  reasonably  close  to  a plane  wave  is  inversely  proportional 
to  the  telescope  aperture  area  (divided  by  r02  ) so  that  there  is  a practical 

limit  of  useful  telescope  size  for  a CENSORING  system.  The  larger  the 
telescope  diameter,  the  longer  we  have  to  wait  for  a good  picture,  with 
waiting  time  increasing  dramatically  as  telescope  size  is  increased. 

For  this  reason,  we  desire  a quantitative  understanding  of  the 
probabilities  involved  in  a CENSORING  system's  operation.  This  paper 
is  aimed  at  the  formulation  of  that  theory.  Here  we  shall  be  concerned 
with  setting  up  the  basic  formulations  and  deriving  equations  suitable  for 
computer  evaluation.  In  a later  paper,  we  shall  carry  out  the  necessary 
computer  calculations. 
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Wavefront  Distortion  Analysis 

The  key  to  analysis  of  the  wavefront  distortion  probabilities  involved 
in  CENSO  iNG  operation  is  the  decomposition  of  a sample  of  the  random 
wavefront  taken  over  the  aperture  into  a set  of  orthonormal  functions  whose 
coefficients  .n  the  decomposition  series  representation  are  independent 
ran  lorn  variables.  From  knowledge  of  the  mean-square  value  of  these 
random  coefficients,  we  can  calculate  the  appropriate  probabilities  for 
CENSORING  system  operation. 

The  orthonormal  decomposition  with  independent  random  coefficients 
is  related  to  the  Ka rhunen - Loe\ e theorem  and  we  can  anticipate  that  the  de- 
velopment of  the  orthonormal  functions  and  the  evaluation  of  the  mean-square 
value  of  the  coefficients  will  depend  on  the  solution  of  a homogeneous  integral 
equation  for  its  eigenfunctions  and  eigenvalues,  respectively.  The  key  to 
that  effort  is  the  development  of  the  kernel  for  the  integral  equation.  The 
kernel  is  developed  from  the  statistics  of  wavefront  distortion  and,  as  we 
m y naturally  expect,  is  related  to  the  phase  struc  ture  function.  However, 
as  we  shall  see,  the  relationship  is  by  no  means  trivial  and  requires  careful 
development. 

In  order  to  calculate  the  kernel,  we  first  have  to  define  in  exact 
terms  the  nature  of  the  wavefront  distortion  statistics  and  the  "portion"  of 
the  distortion  that  is  of  concern  to  us.  We  demote  the  random  wavefront 
distortion  (measured  in  radians  of  phase)  at  a point  r on  the  aperture  plane 
by  0 ( r ) . Over  a circular  aperture  of  diameter  D a random  sample  of 
the  distorted  wavefront  has  a random  average  phase  0 , and  a random 

average  tilt  a . where 

v-  = <4  tt  EFT  1 " dr  W(r , D)  0 (r ) (1 ) 

and 

a = (n  rr  D4)" 1 " d7  W(r , D)  r 0 (r)  , (2) 


. 3 . 


w h c>  r e W{7,  D)  is  an  aperture  func  tion  defining,  in  .he  r-plane,  a circle 
of  diameter  D centered  at  the  origin.  This  aperture  function  is  defined 
by  the  equation 


■ 


II.  if  |7|  < £ D 

w(7.  D)  (3) 

^ 0 , if  | lr  | > £ D 

so  that  in  effect,  it  defines  the  limits  of  integration  for  the  r - integration 
in  Eq.'s  (1)  and  (2),  which  are  otherwise  taken  to  be  over  the  infinite 
r-plane.  The  normalization  in  Eq.  (2)  has  been  chosen  so  that  if  0(r)  = la  • ~r  , 
them  we  would  obtain  from  Eq.  (2)  the  relationship  a = 7 . 

We  note  that  in  forming  a short-exposure  image,  neither  the  average 
phase,  0 , nor  the  wavefront  tilt,  a , disturb  the  resolution  of  the  image. 

The  tilt  produces  an  image  shift  which  is  basically  a nonobservable  in  the 
sense  that  without  special  effort  to  provide  an  absolute  angular  orientation 
reference,  the  effect  of  the  tilt  will  not  be  measurable.  Thus,  from  our 
point  of  view,  the  effective  instantaneous  random  wavefront  distortion  over 
the  aperture  is 

cp(~r;D)  = 0(r)  - 0 - a * r (4) 

If  cp  is  small  enough,  then  the  image  will  be  nearly  diffraction-limited  no 
matter  how  large  0 and  a are.  We  wish  to  calculate  the  probability  dis- 
tribution for  cp  as  the  basis  'or  determining  the  statistics  of  the  operation 
of  a CENSORING  system. 

To  provide  the  basis  for  our  calculations  of  these  statistics,  at 
this  point  we  introduce  the  set  of  functions,  {fa(r;D)]  with  the  orthonormal 
property  that 

( 1 , if  n = n' 

J d7  W(r,  D)  fn*(r;D)  f,,(?;D)  = 

.4  . 


0 , if  n / n 


(5) 


and  the  completeness  property  that  for  any  random  sample  cp(r  D)  we  can 
w rite 


?(r;»)  )3BfB(7;0)  . (6) 

La 

n 

Where  f6„]  is  an  appropriately  chosen  set  of  coefficients.  Because  of  the 
orthonormal  property  of  f , as  defined  by  Eq.  (5),  it  follows  from  Eq.  (6) 
that 

3.  J dr  W(?,  D)  fn*(r;D)  cp(r;L>)  . (7) 

Obviously,  then,  just  as  $ is  a random  function,  3B  is  a random  variable. 
We  also  note  that  since  0 is  a gaussian  random  function,  then  0 and  a , 
being  linear  functions  of  0 are  gaussian  random  variables.  From  this,  in 
turn,  it  follows  that  cp  is  a gaussian  random  function  --  and  this,  in  turn, 
implies  that  3n  , being  a linear  function  of  cp  , is  a gaussian  random 
variable.  This  fact,  together  with  our  ability  to  calculate  the  variance 
of  3,  (which  we  shall  obtain  as  the  eigenvalues  of  the  integral  equation 
defining  fB  ),  will  provide  the  basis  of  calculating  the  probability  of  cp 
taking  a low  enough  wavefront  distortion  form. 

The  key  to  the  definition  of  the  set  of  orthonormal  functions  [f#] 
from  all  possible  sets  of  functions  that  are  orthonormal  over  the  region 
defined  by  W(r,  D)  is  the  requirement  that  the  various  random  coefficients 
3n  must  be  independent.  We  require  that 

( Bnp(D)  , if  n = n ' 

<SB  = < (8) 

jO  , if  n -L  n' 

where  Bn2(D)  denotes  the  variance  of  the  random  variable  Pn  . (In  writing 
B 2 , we  have  chosen,  as  a matter  of  convenience  for  later  work,  to  make  the 

n 

dependence  of  Bb2  on  the  aperture,  diameter,  D , explicit.)  To  see  the 
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implications  of  Eq.  (8),  we  consider  the  quantity 


i = < j*d7'  W(7'.  D)cp*(?';D)cp(?,D)fB(T';D)>  . (9) 

U’e  define  the  covariance  of  cp  as 

Ccp <1  ^ - r ' | ;°)  <cp*(T;D)  w(7';D)>  , (10) 

where  the  homogeneity  and  isotropy  of  the  propagation  statistics  have 
allowed  us  to  write  t h<  dependence  on  r and  7'  in  the  form  of  - "r'| 

We  note  that  by  interchanging  the  order  of  ensemble  averaging  and  inte- 
gration, we  can  rewrite  Eq.  (9)  in  the  for  o 

5 = j’dr'  W(?',D)  (|r  - 7'|  ;D)  f.(r ' ;D)  . (11) 

However,  if  we  use  Eq.  (6)  to  provide  a series  representation  for  cp*(r';D) 
to  be  substituted  into  Eq.  (9),  we  get 

* = <j*  dr'  D)  cp(r;D)  J B.,*  f„V';D)  f„(r';D)>  , (12) 

n' 

which  on  interchanging  the  order  of  integra'ion  and  summation  gives  us 

s = <<p(?;D)~  0,/j*  dr'  W(7',D)fB*(7';D)fB(7';D)>  . (13) 

n' 

The  orthonormal  property  of  f„  as  expressed  in  Eq.  (5)  allows  the  inte- 
gration to  be  performed,  and  the  result  of  this  allows  us  to  reduce  the  sum- 
mation over  n'  to  the  single  term  for  which  n'  = n . Thus  we  get 

6 = <cpfr;D)  gn*)  , (14) 


- 6 - 
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Now  ii  we  again  use  Eq.  (6)  to  provide  a series  representation  for 
cp(r;D)  , we  get 

S = < l > (15) 

n' 

An  interchange  of  the  order  of  summation  and  ensemble  averaging 
leads  to  the  result 

$ = ^ <pn*  0,,)  . (16) 

ft' 

Now  we  make  use  of  the  requirement  that  the  random  coefficients  0n  and 
dn,  be  independent,  as  expressed  by  Eq.  (8),  which  allows  us  to  reduce 
the  summation  on  n'  to  a single  term  with  n'  = n . Thus  we  get  as  a 
consequence  of  the  requirement  of  independence  of  the  P„'s 

S - B,S(D)  f,  (?;D)  , (17) 

which,  when  combined  with  Eq.  (11),  gives  the  basic  Karhunen- Loeve 
homogeneous  integral  equation 

; dr'  W(7',  D)  qp(|r  - r'|;D)  fn(7';D)  - Ba2 (D)  f#(r;D)  . (18) 

The  eigenfunctions  of  this  equation  are  the  orthonormal  functions  we  wish 
to  work  with.  'Ihese  functions  have  statistically  independent  coefficients  in 
a series  representative  of  tp  as  required  by  Eq.  (8).  The  variance  of 
these  coefficients  are  the  corresponding  eigenvalues  of  the  Ka rhunen- Loeve 
integral  equation. 

Our  basic  remaining  tasks  in  this  paper  are  1)  to  see  how  the  prob- 
ability of  the  effective  wavefront  distortion,  cp  , being  adequately  small  can 
be  calculated  from  the  eigenfunctions  and  eigenvalues  defined  by  Eq.  (18), 
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;ind  2)  to  develop  an  expression  for  CT,  to  bo  substituted  into  Eq.  (18). 
However  betore  delving  too  deeply  into  the  se  matters,  we  first  wish  to 
consider  the  question  of  how  we  can  "dedimensionalize"  our  problem  so 
that  a single  case  solution  of  the  Ka rhunen- Loeve  integral  equation  will 
provide  the  basil  treatment  for  all  aperture  diameters  and  strength  of 
turbulence  conditions. 

Uediinension.il  i /.a lion 

The  basic  statistics  of  wavefront  distortion  are  provided  by  the 
wave-structure  function,*  J?(  r)  , where 

r'|  = < |?  (r)  - 0 (r  ')\2  ) . (19) 

It  can  be  shown  that  the  value  of  the  wave -structure  function  may  be 
written  a s 

•Mr)  f>.  88  (r/r0)5/3  , (20) 


where  r0  is  a quantity  with  the  dimensions  of  length.  The  value  of  rQ 
is  determined  by  the  optical  wavelength  in  question  and  the  distribution 
of  the  strength  of  turbulence  along  the  propagation  path.  For  our  purposes 
here,  the  nature  of  that  relationship  is  of  no  consequence.  It  is  sufficient 
to  know  the  value  of  r0  as  this  quantity,  as  vve  shall  see,  completely 
characterizes  the  effective  strength  of  turbulence  for  evaluation  of  the 
performance  of  a CENSORING  system. 

At  this  point,  we  state  without  proof  the  fact  that  we  can  extract 
the  dependence  of  Cjp(|  r - r ' | ; D)  on  both  the  strength  of  turbulence,  as 

* It  should  be  noted  that  although  we  have  spoken  of  0 as  a phase  distor- 
tion, we  actually  consider  it  to  be  a complex  phase  with  its  real  part 
corresponding  to  ordinary  phase  and  the  negative  of  its  imaginary  part 
corresponding  to  log-amplitude  variations.  Thus  both  0 and  the  ran- 
dom coefficients  f3n]  are  complex  quantities.  In  all  of  our  analysis,  we 
have  been  careful  to  introduce  complex  conjugation  where  appropriate, 
though  we  have  not  made  a point  of  the  fai  t that  the  quantities  are  complex. 

In  most  cases,  the  imaginary  part  is  much  smaller  than  the  real  part. 
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defined  by 


. and  the  aperture  diarreter,  D , by  writing 


C0(|r-?'|;D)  = (D/r0f'3  (T(|  x - x'|  ) , (21) 

where 

x r/D  and  x'  r'/D  . (22) 

We  note  in  particular  that  (?(|x  - x'|)  is  independent  of  the  value  of  D . 

The  validity  of  Eq.  (21)  will  be  established  in  a later  section,  where  we 
shall  consider  in  detail  the  evaluation  of  and  will  give  an  explicit 

expression  for  $ . 

If  we  substitute  Eq.  (21)  into  Eq.  (18)  and  change  the  variables 
from  r,  r'  to  x,x'  , noting  that  dr'  goes  into  D2dx'  , we  get  the 
result  that 

; dx'  W(x'.  l)<T(|x  - x'|),3n(x';D)  = B.^DJg.OciD)  , (23) 

where 

3nCx;D)  = f„(Dx;D) 

and 

*n2(D)  1 r2  (D/r0)-B/3B n2(D) 

D-u/3  ro6/u  BnS(D) 

We  note  that  Eq.  (23)  is  a homogeneous  integral  equation  with  eigenfunction 
3n(x;D)  and  eigenvalue  «n2(D)  . However,  since  the  kernel  of  that  integral 
equation,  as  well  as  the  limits  on  the  integration,  is  independent  of  the 
aperture  diameter,  D , then  it  follows  that  the  eigenfunctions  and  eigen- 
values must  be  independent  of  D . We  therefore  may  write  the  eigenfunction 
3,(x;D)  as  r$n(x)  without  any  loss  of  definiteness,  with  the  understanding 


(24) 


(25) 
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that  the  eigenfunction  we  originally  sought,  i.  e.  , fB  (r;  D)  can  be  written 
a s 

f„(r;D)  Pljr/D)  . (2<  ) 

Since  the  eigenvalues  of  Eq.  (23),  i.e.  , {ifln2(D)]  , are  independent  of 
D (and  ol  r0  ),  we  can  without  any  loss  of  information  write  them  as 

^ then  follows  that  the  eigenvalues  we  originally  sought,  namely, 
B„2  ( D)  , c an  be  wr  i tten  as 

B„2(D)  r0_B'3  iPnr  . (2  7) 

We  recall  that  in  accordance  with  the  above  discussion,  and  {f^B  (x)l 

are  the  set  of  eigenvalues  and  eigenfunctions  for  the  Karhunen  - Loeve  homo- 
gen  -ous  integral  equation 

1 dx ' W(x'.  1 ) <F(  j x - >'|  ),3n(x')  2 g„(x)  . (28) 

At  this  point,  we  need  only  develop  an  expression  for  <F(|  x - x'|  ) , in 

accordance  with  Eq.  (21),  to  set  up  the  problem  for  solution  of  the  integral 
equation  in  its  general  form.  Then,  after  obtaining  these  generalized  solu- 
tions, for  any  value  of  aperture  diameter,  D , and  turbulence  parameter, 
rQ  , we  can  obtain  the  eigenvalues  and  eigenfunctions,  fBns(D)}  and 
ff„(r;D)}  from  Eq.'s  (26)  and  (27),  for  that  particular  problem.  In  the 
next  section,  we  take  up  the  evaluation  of  (T  . 

Evaluation  of  the  Kernel 

Though  we  are  interested  in  developing  the  kernel,  ^(|  x - x'j  ) , 

we  shall  proceed  in  that  by  means  of  Ec.  (21),  and  in  most  of  this  section 
shall  be  concerned  with  the  evaluation  of  the  original  kernel,  (]"?-"?'  | ; D) 
as  defined  by  Eq.  (10).  We  shall  find,  at  the  end  of  this  section,  that 
extraction  of  a result  for  S(  | x - x'|  ) then  drops  out  as  a trivial  addi- 
tional manipulation. 


If  we  substitute  Eq.  (4)  into  Eq.  (10),  we  obtain 


r'|;D)=([$(r)-?-a'  T‘"i*[0(r/)'0'n'’  r '] ) 

<0*(7)0(?')>  - (0*r?)0)  - (0f?')0*>  4 <0*0> 

- (0*( 7)  a • 7 ')-<0(7  ')  a*’  7)4(0  * a ' r 'S  + <0  a*  * r) 


. -*  ‘ 
+ < a 


r a 


'> 


(20) 


Making  use  of  Eq.  's  (1)  and  (2),  and  interchanging  the  order  of  ensemble 
averaging  and  integration,  we  can  rewrite  Eq.  (29)  in  the  form 

Qp(|  r - r-|;D)  = <0*(7 )0(7')> 

- (1  tt  D2)-1  “*  d7"  W(7",  D)  <0*(7)0(7")> 

w 

- (In  D2)-1  ’ d 7"  W(7",  D)  (0*(7  ")  0(r  ')) 

* v 

4 (\  rr  D2)" 2 fi*  d7"  d7~  W(7",  D)  W(r  , D) 

X <0*(r  ")  0 (r  ")> 

- (TT  n °4)"  1 " d7~  W(7",  D)(0*(r)  0(7  "))  7"  * 7' 

- ( j1;  n D4)"  1 ” d7"  W(r",  D)(0*{7")  0(7'))  7"  • 7 
4 (f6  n D3)-2  fj  d7"  dr~  Wf?",  D)  W(7"\  D) 

X <0*(7")  0(7"’))  7'  • 7~ 

4 (^g  tt  D3)”2  ""  d7"  dr~  W(7",  D)  W(7~,  D) 

X <0*(7")  0(7")>  r r" 

4 (f4  tt  D4)-SJJ  d7"  d7'"  W(?",  D)  W(7~,  D) 

X <0*(7")  0(7")>(7  • 7")  (7'  • 7"')  . (30) 

We  note  that  in  the  first  four  terms  on  the  right-hand- side  of  Eq.  (30), 
if  the  ensemble  averages  are  each  replaced  by  a constant,  then  the  sum  of 
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the  four  terms  vanishes.  This  means  (hat  we  may  add  or  subtract  a con- 
stant value  from  the  ensemble  averages  in  these  terms.  .Ve  further  note 
that  if  the  ensemble  average  in  any  of  the  remaining  terms  is  replaced  by  a 
constant,  the  integration  is  such  that  the  term  will  vanish.  Hence  we  may 
add  or  subtract  a constant  from  any  of  these  ensemble  averages.  Because 
of  the  stationarity  of  (0*(r)  0(r))  , it  is  identical  in  value  to  (0*(r ')  tf(r’')’1, 

or  (0  (r  )0(r  )),  etc.  We  find  then  that  by  appropriate  manipulation 
(which,  amongst  other  things,  includes  taking  note  of  the  fact  that  for  physical 
reasons  is  real,  so  that  we  can  replace  the  right-hand-side  of  Eq.  (30) 

with  one-half  the  sum  of  its  value  plus  its  complex  coniugate),  we  can  re- 
place each  ensemble  average  of  the  product  of  phases  on  the  right-hand-side 
of  Eq.  (30)  by  minus  one-half  the  ensemble  average  of  the  difference  of  the 


phases  absolute  value  square.  Thus 

<0*(r)0(7')>  a -£^(|r_7'|) 

(31) 

or 

<0*(?~)0(7-)>  = - 5 £ (|7~  - r-|)  . 

(32) 

etc.,  where  £ . the  wave -structure  function,  is  defined  in  Eq.  (19),  with 
the  value  given  in  Eq.  (20). 


now 

rewrite  Eq. 

(30) 

in  the 

form 

M 

r - r'|;D)  = 

T, 

+ Ts  i 

T3  + 

t4  + Ts  + t8  + t9 

. (33) 

T>  = 

- 3 ( |"r  -7 

D 

, 

(34) 

TS  8 

i (i"i 

dr" 

w(7", 

D)  .Ml 

7 -7~\)  , 

(35) 

t3  = 

i (i  n D*)->  P 

u 

d7~ 

W(r-, 

D)  M\ 

7'  - 7~\)  . 

(36) 

t4  = 

i drr  D?)~*J 

T dV 

- d?* 

w(?*. 

D)  W(7~,  D) 

x • 

■ r ** 

D 

t 

(37) 
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x5 

T. 

n 

T„ 


= ? (n  ^tD4)-!  r dr-  W(?-,D)^(|?.  r-|)r'  • r'  , (38) 

£ d7*  W(r-,D)  R(  | r ' - r ' | ) T " • 7 , (39) 

f (£n  dr~  dr~  W(7~,  D)  W(7~,  D) 

X i»(|r"  -7~|)  (7  • 7~)  (7'-  7-)  . (40) 

In  developing  Eq.  (33)  from  Eq.  (29),  we  have  dropped  the  two  terms 

T-  = -$  (71  n D")-2 dr~  d7~  W(7~,  D)  W(7~,  D) 

x *(|7~  -?~|)7'- 7~  , (41) 

Te  * n I)3)"S  IV  d7*  dr*  W(r  *,  D)  W (r~,D) 

1 0 ♦ 

x »(? ~ - r ~|  ) r • 7~  , (42) 

which  correspond  to  the  seventh  and  eighth  terms  in  Eq.  (29).  Our  reason 
for  doing  this  is  that  both  terms  can  be  shown  to  have  zero  value.  To  see 
that  these  terms  do  indeed  vanish,  we  note  that  if  we  consider  pairs  of 
values  of  r and  r that  hold  r “ , r '*  , and  the  angle  between  7~ 
and  r constant,  then  as  we  integrate  about  the  orientation  of  the  7'~-vector, 
everything  in  the  integrand  is  constant  except  r • r *\  As  a result  of  the 
variation  of  this  factor,  the  integration  over  2tt  yields  a zero  value.  (It  is 
interesting  to  note  that  T7  and  T0  correspond  to  the  correlation  of  0 and 
n.  , which,  as  we  should  have  expected,  are  uncorrelated.  ) 

With  Eq.  (33)  established,  we  now  look  into  the  problem  of  using 
Eq.  (20)  for  3 to  simplify  our  expression  for  . It  is  convenient  at 

this  point  to  introduce  the  following  functions: 


&0  (x)  = 3.  44  x6'3 
0,  (x) 


(43) 


iA>  2rr 

3.44  (in)-'  J dx~  x*  j*  dQ  "(x2  + x"s  - 2xx~cos  0>/s  (44) 

0 o'" 

l/a 

8 J dx  ' x-  (x~) 
o 
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2 


(45) 


®3  (X)  = 

IP  2rr 

3.44  (^tt)-1  ^dx'  x'2  n dg  ' cos  9' 
0"  0 

X (x2  + x'a-2xx'  cos  q')P/6  , 

(46) 

* 

II 

1/3 

64  J dx'  x'2  ®a(x') 

(47) 

0 


We  shall  see  shortly  that  can  be  expressed  in  terms  of  these  five 

functions  (or  more  precisely,  three  functions  and  two  constants). 

From  Eq.  (20),  we  see  that 


Tt  = -3.44  (r2+  r'2  - 2rr'  cos  Q' rQ~s/3 

= -3.44  (D  / Iq  J6''3  [(r  / D)2  + (r'  / D)2  - 2(r/D)(r'/D)  cos  9']5/6 
= -(D/r0F3®o(|(?/D)  - (?'/D)|)  . (48) 

Proceeding  in  the  same  vein,  we  see  that 

l/e  2tt 

T2  = 3.44  (in)-»  (D/r0f/3  f dx'  x'  f da' 

o 0 

X [(r/D)2  + x'2-2(r/D)  x”  cos 
= (D/r0^  ®,(r/D)  , (49) 

where  we  have  made  the  transformation  of  variables 

r~  = Dx'  and  dir'  = D2  dx'  , (50) 

ar.d  later  will  also  use 

7'  = Ox'  and  d'r~  = if  dx"  . (51) 
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F or  T3  , working  in  exactly  the  same  way  but  with  r replaced  by 
r'  , we  get 


T3  (D/roF3  ®,  (r'/D) 


(52) 


Working  in  the  same  manner  with  T4  , as  defined  in  Eq.  (3  7), 

we  get 


if?  2tt  1/2 

T4  - -3.44  (i  *r)T^(D/ro)!5/3  " dx*  x"  [*  d9"  n dx~  x~  " d0~ 

0 0 0 0 

X [x  2 + x~*  - 2x"  x'"'  cos  (G"  - 0")]6/s 

1/5  , if?  2n 

- -3.44  ( ^ tt )“ 2 ( D / r0 " dx"  x"  2n  J ? dx~  x~  " d9~ 

o 1 o 

x [x*2+x"2-2x*  x~  cos  ( 0~ ) ]B/6  j.  , (53 

Here  we  have  replaced  the  variable  0~  by  9~+9"  (treating  9"  as  a 
constant  for  the  6~- integral  ion),  and  then  shifted  the  limits  of  that  inte- 
gration from  0 2n49  to  0 « 2~i  . This  then  allowed  the  9"- integration 

to  be  performed,  yielding  a factor  of  2TT  for  the  final  result  in  Eq.  (53). 
We  note  that  the  quantity  in  the  curly  brackets  in  Eq.  (53)  is  directly  re- 
lated to  (x")  . Thus  we  can  write 


1/2 

T4  -8  (D/r0f  '3  1 dx"  x~  @ (x~) 

0 

= - (D/r0/'/3  ®2 


(54) 


For  Ts 
powers  of  ET  1 
into  x • x . 


as  defined  by  Eq.  (38),  we  have  the  necessary  extra 
in  front  of  the  integral  to  convert  r " • r'  in  the  mtegiand 
Following  the  same  procedure  as  above,  this  allows  us  to 
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write 


l 1/2  £TT 

Ts  = 3-44  (7;  nl'MD/iiPJ  dx~  x'  J d0~  x'  x'  cos  (9~- 0') 

X [ (r / D)2  + x~8-  2x  x'  cos  , (55) 

where  the  9'  value  is  referenced  to  the  ~r  orientation  so  that  9*=  0 
when  x'*’  is  parallel  to  r (or  rather  when  ~r~  is  parallel  to  ~r  ).  Q' 
is  similarly  referenced  to  the  r vector.  At  this  point,  we  note  that 
since  we  can  rewrite  cos  (9*-  9')  in  the  form 

cos  (9"-  9')  = cos  9"  cos  9'  + sin  9*  sin  9'  , (56) 

then  we  can  separate  the  9 -integration  in  Eq.  (55)  into  the  sum  of  two 
integrations.  The  first  would  be  multiplied  by  a factor  of  cos  9'  and 
would  have  cos  (9"-  9')  in  the  integrand  of  Eq.  (55)  replaced  by  cos  e*  . 
The  second  would  be  multiplied  by  sin  9'  and  would  have  cos  (9"-  9')  in 
the  integrand  replaced  by  sin  e"  . We  note,  however,  that  since  the  inte- 
grand in  the  second  integration  would  be  odd  in  0"  , the  second  integral 
being  over  the  range  0 *■*  2n  will  have  zero  value.  Thus  we  obtain  the 
result  that 

T8  * 3.44  (nnl*1  (D/r0)P/3  x'  cos  9'  f dx*  x'a  J*  d0*  cos  0' 

o*’  o' 

X [(r/D)®  + x"a  - 2x  x"  cos  e']**5 

* (D/r0P*  (r'/D)  cos  9'  0,  (r/D)  . (57) 

Exactly  the  same  procedure  applied  to  the  evaluation  of  T,  leads  to 
the  result 

Ta  * (D/r0f/3  (r/D)  cos  9' 0,  (r'/D)  . (58) 
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evaluation  of  Tg  follows  the  same  general  approach  as  the 

above.  In  this  case,  we  have  enough  extra  powers  of  Dr  i in  front  of  the 
integral  to  allow  us  to  convert  (?  - ?')  (?'  • ?-)  into  [(7/D)  - ?][(?'/!)). 


Thus  we  can  write 


2tt 


2tt 


TS  -3-  44  (S4”)“2  J dx  X*  " de"  ’ dx*  X*  d0~'  (r/D)(r'/D)  x~  x 

° 0 o"  o 


cos  9-  COS  (9~  - 9')  [ x~s  4 x~2  - 2x"  X~  COS  (e'~-  0-)]B/6 


= -3.44  (C  nr*(r/D><r*/D)  fV  x„a  Vfx_  ^ ^ 


X e»  [8-4(3*.  9')]l»'s4»*i.2«-  x-  cos  t~fn 


(59) 


In-order  to  obtain  the  final  form  of  Eq.  (59).  we  have  replaced  e~  by 

0"+  (trea“ng  38  3 C°nStant  for  th«  9*"- integration),  and  then  adjusted 
tbe  limits  of  the  new  9~-integration  from  9"  ~ 2n  + 9~  to  0 „ 2n  . Now  we 
can  rewrite  cos  [9~  + (g—_  e')~  as 


9 + (9  - 9 )]  = cos  9~  cos  (9'- 9')  + sin  9~  sin  (9'- 0')  , (60) 


and  using  the  same  arguments  as  were  used  to  develop  Eq.  (5?)  from  Eq.  (55), 

^droPPmg  thc  sin  9"  dependence  because  it  leads  to  an  odd  integrand 
>n  Q*'  , we  obtain 


T9  = -3.  44  (^  -r)-a(r/D)(r  VD)  'jx-x-a 

o' 

2n 


2tt 


, l/s 

X | ‘ dx~  -"2 


X * u d9"'  cos  9"  (x'2  + x“’2-2x'  x'*' 
0 


J*  d9~  cos  0*  cos  (9".  e') 

9~F/6  }.  (61) 


X cos 


First  of  all,  we  note  the  close  relationship  between  th 


e quantity  in  the  curly 


brackets  in  Eq.  (61)  to  the  integration  defining  ®3  in  Eq.  (46).  Second, 


we 
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note  that  the  only  ©"-dependence  in  Eq.  (61)  is  cos  Q~  cos  (9~  - Q' ) 
Since 

2n 

J de~  c°s  Q~  cos  (e~-  0')  = n cos  0' 


(62) 


we  see  that  Eq.  (61)  can  be  rewritten  as 


T9  = - (D/rg)6^  (r/D)(r'/D)  cos  0' ffl4  . (63) 

With  all  of  these  results  in  hand,  we  can  now  rewrite  as 

given  by  Eq.  (33)  in  the  form 

Cv(l  r-r'|  ;D)  = (^(r  , r',  Q ; D) 

= (D/re^(-^(|(7/D)-(?7D)|)  + (r/D) 

+ Oj  (r'/D)  - Qfe  + (r'/D)  cos  6'  ®3(r/D) 

+ (r/D)  cos  0'  Qfe  (r'/D)  - (r/D)(r'/D)  cos  0'  ©,  } . (64) 

In  writing  Eq.  (64),  we  have  taken  the  liberty  of  introducing  the  notation 
C<P(r  , r',  9'  ; D)  to  make  it  explicit  that  the  dependence  on  |r-r'|  can 
just  as  well  be  considered  a dependence  on  r,r'  , and  0'  . We  note  that 
although  Q'  has  been  defined  in  terms  of  the  7'  , physically,  and  in  our 

future  work,  it  can  be  considered  as  simply  being  the  angle  between  7 and 
r ' , if  we  choose. 

The  first  thing  we  wish  to  note  in  considering  Eq.  (64)  is  the  fact 
that  it  is  indeed  a function  of  7/D  and  r'/D  multiplied  by  ( D/r 
thus  justifying  the  assertion  that  led  us  to  write  down  Eq.  (21).  In  fact, 

comparing  Eq.  's  (21)  and  (64),  we  see  that  the  value  of  ®(|x-x'|)  can 
be  written  as 
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<S(|x-x'|)  - *(x  , x',  9') 

- ®0  (|  x-x'l  ) + @, (x) + 0, (x')  - ©g 

+ x'  cos  9'  (J^  (x)  + x cos  0'  @3  (x')  - x x'  cos  0'  @4  . (65) 

With  this  expression  in  hand  for  K , with  the  ©-functions  defined  by 
Eq . ' s (4  3)  to  (47 ),  we  can  now  turn  our  attention  back  to  the  problem  of 
solving  the  integral  equation  for  the  eigenvalues  and  the  eigenfunctions 

3,  (x)  , in  accordance  with  Eq.  (28). 

Integral  Equation  Reduction 

.As  indicated  at  the  start  of  this  paper,  we  shall  not  attempt  here  to 
solve  the  Karhunen- Loeve  homogeneous  integral  equation,  the  pertinent 
form  of  which  is  given  by  Eq.  (28).  Ultimately  we  intend  to  solve  this  equa- 
tion using  numerical  techniques,  in  particular  the  Givens-Householder 
method.  At  this  point,  however,  we  are  interested  in  numerical  procedures 
that  will  simplify  this  equation.  We  note  that  the  basic  integral  equation  is 
two-dimensional  and  that  as  a i onsequencc,  the  size  of  the  matrix  required 
to  obtain  any  reasonable  resolution  over  the  aperture  for  our  eigenfunctions 
will  be  unreasonably  large.  We  propose  to  avoid  this  problem  by  introducing 
a separation  of  variables  in  the  eigenfunction. 

We  postulate  that  the  eigenfunction  (x)  can  be  separated  into  a 
radial  dependence  and  an  azimuthal  dependence,  and  further  postulate  that 
the  azimuthal  dependence  can  be  written  in  the  form  exp  (i  q 9)  , where 

q-0,±l,i2,.,,  is  a quantum  number'  for  a set  of  solutions  cor- 
responding to  a subset  of  the  n ' quantum  numbers.  " The  radial  depen- 
dence would  have  its  own  set  of  "quantum  numbers",  p = 1 , 2 , 3 

with  each  combination  (p , q ) corresponding  to  an  element  in  the  set  that 
was  "counted"  by  n . Thus  we  would  write 

(x)  = S’(x)  exp  (i  q 9)  , (66) 
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where  the  superscript  q over  the  radial  dependence  term,  i.  e.  , the 
function  ft  , is  used  to  indicate  that  for  each  value  of  q , we  may  ex- 
pect to  obtain  a different  set  of  radial  dependence  functions.  The  eigen- 
value would  now  be  written  as  ^ in  place  of  SP#C  . 

To  validate  our  hypothesis  concerning  the  validity  of  Eq.  (66),  we 
shall  substitute  Eq.  (66)  for  (x)  into  the  Karhunen-Loeve  homogeneous 
integral  equation  given  by  Eq.  (28),  and  show  that  it  leads  to  self-consis- 
tency in  the  sense  that  the  form  of  ^ (x)  will  be  found  to  have  .the  form 
given  by  Eq.  (66).  While  this  is  not  a rigorous  proof  of  the  validity  of  the 
separation  of  variables,  our  manipulations  will  define  the  key  steps  re- 
quired to  develop  such  a rigorous  proof.  We  shall  not  concern  ourselves 
further  with  the  matter  of  a rigorous  proof. 

If  we  substitute  Eq.  (66)  for  HjB(x)  into  Eq.  (28),  we  obtain 

dx'  W(x',  1 ) d(|  x-x'  |)  gn  (x') 

1/2 

' dx'  x'  r de'  S(x.x',e'-e)RMx')  exp  (i  q 0')  , (67) 

0 0 

where  here  we  have  chosen  an  arbitrary  angular  reference  point  so  that 
we  can  define  angles  0 and  9'  associated  with  x and  x'  , respectively, 
rather  than  merely  having  9'  defined  as  the  angle  between  x and  x'  . 
Then  in  writing  d , we  took  note  of  the  fact  that  the  ©'-dependence  indi- 
cated in  Eq.  (65)  was  in  this  case  a dependence  on  9'-  9 . Now  if  we 
replace  9'  in  Eq.  (67)  with  9'f  9 and  then  readjust  the  limits  of  the 
9'- integration  from  9*-.  2-  + 9 to  0 ^ 2n  , we  obtain 

f dx'  W(x',  1)  d(|  x-x'| ) 3n(x') 

1/2 

= exp  (i  q 9)  "dx'  (x,  x' ) Rj*  (x')  , (68) 

0 
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where  we  have  used  ft  to  mean 

2r, 

R,(x,x')  = x 'do'  s(x,  x',9')  exp  (i  q 9')  . (69) 

o' 

Now  combining  Eq.  's  (68)  and  (28),  we  see  that 

. 1/2  % 

3„  G)  = exp  (i  q 9)  ^-2  J dx'fi,(x,x')SpMx')^  . (70) 

o 

Comparison  of  Eq.  (70)  with  Eq.  (66)  shows  that  our  assumption  of  Eq.  (66) 
for  ,jb(x')  leads  to  se  If- c onsistent  re  suit  s for  gn  (x)  • Combining  Eq.  1 s 
(66)  and  (70),  we  obtain  the  equation 

1/2 

J dx' ft,  (x,  x')  S’ (x')  = ¥p2,  Sp’(x)  . (71) 

o’ 

This  is  a Karhunen- Loe ve  integral  equation  type  definition  for  the  radial 
function,  S’  , with  a kernel,  ft,  , which  can  be  different  for  each  value 
of  the  'quantum  number,"  q . 

Eq.  (71)  provides  us  with  a basis  for  calculation  of  the  complete 
set  of  eigenvalues,  , (or  SBb2  in  our  original  notation),  and  together 
with  Eq.  (66)  provides  a definition  of  our  two-dimensional  eigenfunction, 

S’  (x)  exp  (i  q 0)  (or  f),  (x)  in  our  original  notation).  The  set  of  kernels, 
{ft,}  , can  be  obtained  by  substituting  Eq.  (65)  into  Eq.  (69).  We  get 

2n 

ft  0 (x,  x' ) = -x  ” d9'  ®0  ([  x2  + x'2  - 2x  x'  cos  9']1/2 
0 

+ 2n  x'  [«,  (x)  + ®,  (x')  - «8]  , (72) 

2tt 

ft±  , (x,  x' ) = -x'  ',d0'®o([x2  + x'2  - 2x  x'  cos  0']1/2)  cos  (00 

0 

+ TT  x'  [ x'  ©3  (x)  + x ®3  (xO  - X x'  ®4  ] , (73) 
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and  for  the  magnitude  of  q greater  than  one 

2n 

ft?(x,x')  = - x*  ’ dO'  @0  ([  x2  + x2  - 2x  x cos  8']1/2)  cos  (q  ') 

0 

for  q ±2  , ±3  , ±4  , . . . (7 }) 

In  obtaining  these  results,  we  have  made  use  of  the  fact  that  exp  (i  q 9') 
can  be  written  as  cos  9'+  i sin  9'  , and  noted  that  the  sin  0'  leads  to 
integrands  odd  in  £>'  so  that  their  value  after  integration  of  0 2n  vanishes. 

The  combination  of  Eq.  's  (26),  (27),  (43),  (44),  (45),  (46),  (47),  (66), 
(7  1 ),  (72),  (73),  and  (74)  provides  the  basis  for  calculating  the  eigenvalues 
and  eigenfunctions  we  shall  need  to  determine  the  probability  of  an  accidental 
occurrence  of  a low  energy  wavefront  distortion  condition  so  that  the 
CENSORING  system  will  be  able  to  produce  a near  diffraction-limited  image. 
In  the  next  section,  we  take  up  the  problem  of  calculating  this  probability, 
given  the  set  of  eigenvalues.  As  noted  before,  we  leave  the  problem  of 
numerically  evaluating  the  eigenvalues  and  eigenfunctions  for  treatment 
in  a subsequent  paper. 

Probability  Formulation 

The  key  to  the  evaluation  of  the  probabilities  associated  with  a 
CENSORING  system's  performance  is  to  recognize  that  this  is  essentially 
equivalent  to  a study  of  the  probability  distribution  of  the  effective  mean- 
square  wavefront  distortion  over  the  aperture.  The  term  "effective"  as 
used  here  refers  to  the  fact  that  we  are  only  interested  in  wavefront  varia- 
tions excluding  tilt  and  average  phase  variations  --  i.e.  , the  effective 
wavefront  distortion  is  to  be  calculated  from  cp(r;D)  and  not  from  0(r)  . 

We  write  the  mean-square  wavefront  distortion  over  the  aperture  as 

£2=(iTTD2ri  !*  dr  W(r;D)  |<P(r;D)  |2 


- 22  - 


(7b; 


It  should  bo  recognized  that  the  term  "mean"  in  reference  to  A2  refers 
to  an  average  over  the  aperture  and  not  to  an  ensemble  average.  Thus 
just  as  the  effective  wavefront  distortion,  cp  , is  a random  function,  we 
see  from  Eq.  (75)  that  A2  is  a random  variable.  If,  at  some  instant, 

A1  is  small  enough,  then  we  may  expect  nearly  diffraction-limited  quality 
for  an  image  formed  at  that  instant.  The  problem  we  face  in  calculating 
CENSORING  system  performance  is  one  of  calculating  the  probability  that 
A“  will  be  sm;ll  enough.  We  set  as  a nominal  threshold  the  requirement 
that  A"  < r a d ia  n - squa  r e as  the  dividing  line  between  good  and  poor 

images.  Our  problem  is  to  calculate  the  probability  of  A2  being  less 
than  A T 2 , this  being  the  probability  that  the  CENSORING  system  will 

see  good  enough  conditions  to  allow  an  image  to  be  formed. 

If  we  substitute  Eq.  (6)  into  Eq.  (75),  we  get 

A2  =(ln  l J*  d?  W(?,D)f#V)fn,(?)  . (76) 

n.  »' 

Now  making  use  of  the  orthonormality  of  fn  , as  defined  in  Eq.  (5),  we 
can  reduce  Eq.  (76)  to  the  form 

a2  = (|  n o?r‘  ^ e„*  a, 

D 

= l [B„/{iTT  D2)]  . (77) 

a 

Thus  the  mean-square  wavefront  distortion  is  seen  to  be  the  sum  of  the 
square  of  a set  of  gaussian  random  variables,  B„/(x  n D2)1/s  . We  recall 

that  according  to  Eq.  (8),  the  random  variable  has  a variance  given 
by  the  eigenvalue  Bn2(D)  , so  that  the  variance  cf  3b/(t  n D2)l/2  can  be 
written,  using  Eq.  (27),  in  the  form 
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Var  [0,/(J  n I?)1'2]  = (intf)-i  B„2(D) 


= ( D / r0 p 3 (4/tt)  flB2  . (78) 

For  convenience,  we  shall  denote  this  variance  by  a 2 (or  a.  2 ) as 
appropriate. 


aB2  = Var  [0B/(|  tt  I^^2  ] . (79) 

If  we  can  compute  the  eigenvalues,  iPn2  (or  iflp  2 ) for  the  dedimension- 

alized  Karhunen - Loeve  Homogeneous  Integral  equation  of  Eq.  (28)  [ or  of 
Eq.  (71)],  then  we  can  immediately  write 

CB2  = (D/r0f/3  (4/TT)!flB=  , (80) 

or 

Jr.,2  r (D/r0^  (4/n)3V,2  • (80') 

It  now  follows  that  the  probability  of  CENSORING  system  at  any 
instant  seeing  low  enough  distortion  to  allow  an  image  to  be  formed  can 
be  written  as 


Pc.B|or  = Prob  (CENSORING  System  Forming  an  Image) 

= Prob  (A2  s At2)  . (81) 

Since  the  random  variables  0#  are  independent  and  gaussian  distributed 
with  variance  aB2  , it  follows  that 

ao 

Pc.B.or  = ^ (2n  a.s)“1/s  ’ d\  exP  (-i  x.a/°»2)  - (82) 

»=  I J 

Lt.lt. 

where  the  limits  on  the  integration  are  to  be  understood  as  a composite 
limit  on  the  product,  or  rather  on  the  n-tuple  multiple  integral.  The  limit 
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corresponds  to 


Limit 


( L **'  s 4'= ) 


(83) 


Eq.'s  (82)  and  (83),  in  concept  .at  least,  provide  a basis  for  the  calculation 
o{  Pc,nl9r  • >*e.,  the  probability  of  a CENSORING  system  image  being 
formed.  However,  because  of  the  infinite  limits  on  n in  these  two  equa- 
tions, no  practical  calculations  can  be  performed. 


To  provide  a practical  basis  for  carrying  out  the  calculation  of 
Pc  ■ n ( or  > vve  need  to  truncate  the  series.  To  do  this,  we  first  note  that 
in  accordance  with  results  which  we  have  previously  obtained  elsewhere,5 
we  know  that  the  ensemble  average  value  of  A2  can  be  written  as 


( A2)  = 0.  1345  ( D/r0f/3  . (84) 

Now  if  we  assume  that  the  eigenfunctions  are  arranged  in  such  an  order 
that  cB2  is  monotonically  decreasing  with  n , then  we  know  that  if  N 
is  large  enough  so  that 

N 

l c.2  <A2>  - e ATb  . (85) 

n ~ I 

then  if  e has  been  chosen  to  be  a small  enough  quantity,  we  may  con- 
sider N , rather  than  ® , to  be  the  practical  upper  limit  on  the  n 
dependencies  in  Eq.  's  (82)  and  (83).  As  a practical  matter,  we  would 
replace  A T 2 in  Eq.  (83)  with  AT2(l-e)  . 

This  has  the  effect  of  saying  that  above  some  value  of  n (namely 
n = N ),  we  are  not  particularly  concerned  with  the  exact  amount  of  wave- 
front  distortion  introduced  by  each  degree  of  freedom.  The  exact  value 
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of  those  9,'s  does  not  concern  us  since  we  know  that  the  corresponding 
variances,  sB2  , are  so  small  that  the  values  of  the  9„'s  will  be  toler- 
ably small.  We  expect  the  contribution  of  all  those  higher  order  terms, 
which  we  are  suppressing,  to  the  mean-square  wavefront  distortion  to 
only  be  of  the  order  of  eAT2  , which,  by  our  choice  of  e . we  have  made 
sure  is  tolerably  small. 

In  an  actual  calculation  of  Pc<#  r , we  would  fiist  compute  an 
ordered  series  of  eigenvalues  $B2  . Then  using  Eq.  (80),  we  would  com- 
pute the  variances,  ~B2  . By  applying  Eq.  (85),  we  could  then  determine 

the  truncation  level,  N . At  that  point,  our  calculation  would  reduce  to 
the  evaluation  of  the  integral 


where 


i . . • d x.  "i  \ 


; dx,  dxg.  . .dx,  1 1 (2H  ~b2)-i/2  exp  (-i  xB2/oB2)  . (86) 

n=  1 


.imit  s ( ^ x.2  s A T2 (1  - e)  j 


Our  problem  is  thus  reduced  to  the  numerical  evaluation  of  the  integral 
in  Eq.  (86),  having  first  solved  the  Karhunen- Loeve  homogeneous  integral 
equation  for  the  eigenvalues.  Neither  is  a trivial  numerical  task,  but  in 
practice  can  be  expected  to  be  rather  straightforward,  if  somewhat  pon- 
derous in  terms  of  required  computer  effort.  These  tasks  are  taken  up  in 
Part  II  of  this  report. 
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PART  II 


Numerical  Evaluation  of  Pro 
Governing  the  Performance 
CENSORING  System 


abilities 
of  a 


2 


Introduction 

In  Pc'  rt  I of  this  report,  a formal  basis  was  developed  for  the  analy- 
sis of  the  expected  performance  of  a CENSORING  system.  The  basic  quantity 
of  interest  was  identified  as  the  probability  that  at  any  instant  of  time,  the 
wavefront  distortion  over  the  entrance  aperture  would  have  an  rms  devia- 
tion from  a plane  of  less  than  one  radian.  (The  term  " rms  " is  used  here 
in  the  sense  of  an  average  over  the  aperture.  ) Because  the  CENSORING 
system  is  designed  to  form  a short  exposure  image,  the  wavefront  deviation 
is  to  be  measured  relative  to  the  optimally  chosen  tilted  plane,  i.e.  , that 
plane  whose  tilt  is  such  as  to  minimize  the  rms  deviation. 

It  was  shown  that  the  probability  of  interest  could  be  calculated 
in  terms  of  a multi-dimensional  gaussian  distribution  in  a hyper-space, 
where  the  hyper-space  was  defined  in  terms  of  a set  of  functions  which 
could  be  used  to  decompose  a sample  of  the  randomly  distorted  wavefront 
taken  over  the  aperture  into  a set  of  statistically  independent  components. 
Because  the  wavefront  distortion  is  a gaussian  random  process,1  and  be- 
cause the  magnitude  of  each  of  these  independent  random  components  is 
obtained  by  a linear  process  from  the  random  wavefront  distortion,  it 
follows  that  each  component  is  a gaussian  random  variable.  The  random 
amplitude  of  each  component  represents  one  of  the  dimensions  in  the  hyper- 
space, and  the  probability  of  interest  is  the  probability  that  all  of  the  ran- 
dom variables  will  take  on  small  enough  values  at  some  instant  of  time. 

Because  the  functions  used  for  the  decomposition  of  the  wavefront 
are  a set  of  functions  that  are  orthonormal  over  the  space  defined  by  the 
CENSORING  system's  aperture,  it  follows  that  the  mean  square  deviation 
of  the  wavefront  at  any  instant  (with  the  mean  taken  as  an  average  over  the 
aperture)  is  just  equal  to  the  sum  of  the  square  of  the  component  amplitudes 
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in  thi1  wavi'front  decomposition.  This  means  tliat  the  probability  of  the 
rms  wavefront  distortion  being  less  than  or  equal  to  one-radian  is  just 
the  probability  that  the  random  components  will  define  a point  in  the  hyper- 
space  that  lies  within  a hyper-sphere  of  one -radian  radius  and  centered  at 
the  origin.  Since  the  components  each  obey  an  independent  gaussian  distri- 
bution, the  problem  of  evaluating  the  probability  of  interest  can  be  seen  m 
reduce  to  a multi-dimens  ional  integral  within  a unit  hyper-sphere  of  a set 
of  gaussian  distributions.  All  we  need  in  order  to  be  able  to  carry  out  this 
evaluation  is  information  on  the  variances  to  be  associated  with  each  element 
of  the  set  of  gaussian  distribution.  In  the  previous  work,  it  was  shown  that 
these  variances  c ould  be  obtained  by  solving  the  Karhunen - Loeve  integral 
equation  associated  with  the  wavefront  distortion  statistics,  and  that  they 
were  proportional  to  the  eigenvalues  of  that  equation. 

The  single  wavefront  distortion  statistic  required  for  this  problem 
is  the  wave  structure  function,  J?(r)  , where 

Mr)  6.  88  (r/r0F3  . (1) 

The  orthonormal  function  set  ffn("r)]  and  the  associated  set  of  variances, 
fct/l  , which  appropriately  decompose  the  distorted  wavefront  and  which 
define  the  gaussian  probability  distributions  of  interest  in  our  hyper-space 
integration  have  been  shown  to  correspond  to  the  eigenfunctions  and  to  be 
proportional  to  the  eigenvalues  of  a Karhunen- Loeve  integral  equation  util- 
ising a function  of  ft(r)  as  the  kernel  of  the  integral.  The  relationship  be- 
tween the  set  of  variances  of  interest  and  the  set  of  eigenvalues,  f Bn2 } , is 

giver,  by  the  equation 

* B„=/<  1 „ D’)  . ,2) 

In  seeking  a solution  of  this  integral  equation,  it  has  been  shown 
that  the  set  ffn(r)]  can  be  decomposed  in  subsets  by  separation  of  variables 
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into  polar  coordinates,  i.  e.  , 


r = (r,9) 


(3) 


It  has  been  shown  that  we  can  write 

fn(7)  - Rp’fr)  exp  (1  q 6)  , q = -«,  ...-2,  -1,0,  1,2,...+® 

p = 1 , 2,  3,  . . . +®  , (4) 

where  the  function  R^r)  represents  a set,  on  p , of  functions  that 
satisfy  a homogeneous  integral  equation  with  a kernel  that  is  different  for 
each  value  of  q . (Actually  the  kernel  for  q and  -q  are  identical.)  The 
eigenvalues  of  this  integral  equation,  B?  (j3  can  be  equated  with  the  eigen- 
values of  the  original  integral  equation,  i.e.  , 


3 = B 2 

P»  <1  n 

and  equivalently,  the  variance  of  interest  can  be  written  in  the  form 
where 


(5) 


p r i 


°p. 


2 

q 


(6) 


The  integral  equation  defining  R^  (r)  and  1^,  * involves  the 

aperture  diameter  of  the  CENSORING  system,  D , and  the  basic  wavefront 
distortion  turbulence  parameter,  r0  . It  is  convenient  to  cast  the  integral 
equation  in  a form  which  is  independent  of  these  two  parameters  so  that  a 
single  set  of  numerical  solutions  to  the  integral  equation  can  be  applied  for 
all  possible  values  of  D and  rQ  . It  has  been  shown  that  if  we  define  the 
eigenfunctions,  3?^  (x)  and  q2  by  the  integral  equation 

1/2 

f dx'fiq(x,x')  = 93p<2VM  ■ (7) 

o 
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then 


V -= 


DU/3  r -6/3  m 2 


and 


f.(r)  = V<r)  exP  (i  9 6)  = Hp«(r/D)  exp  (i  q 9) 


(8) 


(9) 


2tt 


If  we  define  qjx.x')  as 

(x,  x ) = - X d9  ' gfc  (r  X2  4 x'2  - 2xx'  cos  e ']1/2)  cos  (q0  ') 
0 

then  the  kernel  for  the  integral  equation  of  Eq.  (7)  can  be  written  as 


«,(x.x')  =«q(x.x')  if  q = ±2,  ±3,  ±4 


XX 


<M*.  O = «,(x.  x')  4 T7  X '[  x'&3(x)  4 X03(x')- 
q0(x,x')  = ft0(x,x')  4 2nx'[0j(x)  tQ^x')  - Q ] 


(ID 

(12) 

(13) 


The  0 -functions  are  defined  as 


% (x)  3.  44  x£/3 


1/2  2tt 

@,(x)  - 3.44  (i nr1  Jdx'  x'  ’ de"(xz4x'2-2xx' 

o'*  o' 

1/2 


i)2  = 8 [ dx"  x~  0,(x')  , 


(14) 

cos  0'*')6/6, 

(16) 


1/2 


2tt 


®3  (x)  - 3.44  (in)"1  * dx"  x"2  " d0  * cos  0'(xs4x'2 
o"  o'- 

- 2x  x*  cos  0"  J6'6 
1/2 

= 64  r dx ~ x~2  03(x") 

0 


(17) 

(18) 
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. (10) 


(15) 


— 


We  note  in  passing  that  if  we  had  been  interested  in  a system  in 
which,  even  during  a short  exposure,  wavefront  tilt  had  to  be  considered 
a portion  of  the  wavefront  distortion  (which  it  does  not  have  to  be  lor 
normal  short  exposure  purposes),  then  the  only  change  in  the  above  results 
would  have  been  to  modify  Eq.  (12)  to  the  form 


q^x,  x')  = tf,(x,  x') 


(12') 


Our  basic  problem  is  to  solve  the  integral  Eq.  (7)  for  all  of  the 
eigenvalues,  (3^  * } and  then  for  a particular  set  of  values  of  D and  r0 
obtain  the  corresponding  set  of  variances  { 52}  in  accordance  with 
Eq.  's  (2)  and  (8),  from  which  we  get 


p.  q 


4 /_D  f/3 
tt  V r o ) 


(19) 


The  second  half  of  the  problem  is  the  evaluation  of  the  integral  in  hyper- 
space which  defines  the  probability  that  at  any  instant  the  random  wave- 
front  distortion  relative  to  the  optimally  chosen  tilted  plane  will  have  a 
mean  square  value,  averaged  over  the  aperture,  of  one  radian  squared 
or  less.  This  probability  can  be  written  as 

PctNSO*  (2TT  J dXP.,  eXP  ('a  ’ (20) 


ft  q 


sph 


where  the  limits  on  the  integral  correspond  to  a hyper-sphere  for  which 


P#  4 


xp. 


(21) 


In  Eq.  's  (20)  and  (21),  the  product  and  the  summation  run  over  all  possmle 
combinations  of  values  of  p and  q . 


1 

I 
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Iii  the  next  sec  tion,  we  take  up  the  problem  of  casting  the  integral 
equation  ot  Lq.  (?)  in  a form  suitable  for  numerical  evaluation.  In  the 
sections  after  that,  we  shall  first  present  the  numerical  solution  technique 
and  results,  and  then  go  into  the  problem  of  formulating  an  explicit  form 
for  the  hyper-space  integral  of  Eq.  (20)  for  various  values  of  D/rn  . Then 
we  shall  move  on  to  consider  the  evaluation  of  that  hyper-space  probability 
integral.  Based  on  the  results  of  this  evaluation,  we  shall  present  a gen- 
eral discussion  of  the  expected  performance  of  a CENSORING  system. 

Integral  Equation  Numeric  al  Formulation 

Ihe  numerical  solution  of  the  Karhunen-Loeve  integral  equation 
presented  in  Eq.  (,)  can  be  developed  using  standard  numerical  techniques. 
However,  it  will  greatly  simplify  our  numerical  treatment  if  we  first  recast 
that  equation  into  an  equivalent  form  in  which  the  kernel,  i.e.  , ft  (x,x')  is 
replaced  by  a kernel  that  manifests  symmetry  between  x and  x'  . [We 
note  that  the  leading  factor  of  x'  in  the  right  hand  side  of  Eq.  's  (10),  (12), 
and  (13)  dest  r oys  the  symmetry  of  ft  ^(x,  x' ) . ] 

In  order  to  obtain  the  desired  symmetry,  we  introduce  the  following 
symmetrizing  functions: 

2n 

fl,5  (x,x')  = -(xx'F2  ’ d9'(J)0([x2  + x'2-2xx'  cos  0']^)  cos  (qB')  , (22) 


o 

ft0'J  (x,  x')  = ft0s  (x,  x')  + 2rr  (x  x')1/2[ffl,(x)  + @,(x')  - ©2  ] , (23) 

S,s(x,x')  =ft,S(x,x')  + TT  (xx'p^lx)  +«3(x')  - ffl*  ] , (24) 

ft_!?  (x,  x' ) = ft  J5  (x,  x' ) , (25) 

5V(x)  = x1/2Sp’W  . (26) 


Now  it  follows  from  direct  substitution  that  Eq.  (7)  can  be  rewritten  as 
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(27) 


’ dx'  «,5(x.x')  = SPPl,2  *SpMx) 

0 

which  integral  equation  lias  a symmetric  kernel  and  can  be  solved  using 
straightforward  numerical  techniques.  The  eigenvalues  of  Eq.  (27)  are 
identical  to  those  of  Eq.  (7),  and  the  e igenfunc  t ions  of  Eq.  (7),  i.  c.  , \q  (x) 
can  be  obtained  from  the  c*  igenfunc  tions  of  Eq.  (27),  i.e.  , * 9lpq(x)  by  use 

of  Eq.  (26). 

To  obtain  the  eigenvalues  and  eigenfunctions  of  Eq.  (27),  it  is  nec- 
essary to  transform  the  integration  into  a summation,  thereby  obtaining  a 
homogeneous  set  of  simultaneous  equations  with  a determinant  that  deter- 
mines the  eigenvalues  and  eigenfunctions.  In  order  to  replace  the  integra- 
tion by  a summation,  we  subdivide  the  range  x = 0 to  0.  5 into  20  sections 
and  consider  the  values  of  x at  the  midpoint  of  each  section,  which  we 
denote  by  Xj  for  i = 1 to  20  . We  can  then  make  the  replacement 


f/2dx'  (x , x * ) sV(x') 

o 


2 O 

t . 

1=  l 


Kqs(i.i')  SV(‘') 


(28) 


which  allows  us  to  rewrite  Eq.  (27)  as 
20 

Kq  s ( i . i ) 5 V(i')  = *P.q2  5 V(i) 

l '=  l 


where 

jy(i)  = 5v(xi) 

and 

KqMi.i')  = (i/20)  ttqMXj.X,,) 


(29) 


(30) 

(31) 


Eq.  (29)  presents  us  with  the  straightforward  problem  of  obtaining 
the  eigenvalues  and  eigenfunctions  of  the  matrix  Kq5(i.i')  . This  is  a 
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standard  problem  in  numerical  analysis.  The  20  x 20  size  of  the  matrix 
means  that  we  are  dealing  with  a rather  small  problem  as  such  computa- 
tions go  on  a large  high-speed  digital  computer.  The  mathematical  pro- 
cedure that  we  shall  use  is  based  on  the  well-known  Givens -Householder 
algorithm®  and  associated  procedures.  (The  details  of  the  actual  compu- 
tation will  not  be  discussed,  as  the  computer  program  utilized  a proprietary 
CDC  subroutine  for  this  task.  ) 

Computationally,  the  evaluation  of  Kq*  (i,  i'J  matrix  elements 
represented  almost  as  large  a task  as  the  determination  of  the  eigenvalues 
and  eigenfunctions  of  this  matrix.  The  pertinent  expressions  can  be  written 
as 


JV(i.  »')  = 6.88  (.  02 5)(Xj 


1/2 

+ 8 f du  u P d 9 


' d®  (xi2+  x,,a-  2x,  x,,  cos  9FS 
*•  d 

[ (x,8  + u2  - 2xt  u cos  9)^/6+  (xs  8 + u2 


0 0 


- 2 Xj  ,u  cos  9f/8  ] 


1/2  4/2 

_ 64  du  u dv  v dO  (u2  + v2  - 2uv  cos  Q)6^6 

* * w 

0 0 0 


(32) 


Kt,s(i.i')  = 6.88  (.025)(x,  x, , F 


2 - d9  (x,2  + x,  ? - 2x,  x,  , cos  9Fscos  9 
o 


1/2 


+ 64  du  u3  ' d0  [xt(x,  8+ua-2x,,u  cos  ef/6+  x,,(x,8 


u2  - 2 Xj  u cos  9 F8  ] cos  9 


1/2 


1/2 


- 4096  x,  x,  , du  u2  ' dv  v2  ' d9  (u2  + v8  - 2 u v cos  eF^cos  0 v , (33) 


and 


K,*  (i.  i')  = -6.  88  (.  02 5)(x,  x,,Fs  ' d0(xt8  + 


2x,xt,  cos  9f/6cos  (q9 ) , 
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for  | q | > 1 


(34) 


Wt*  have  utilized  the  trapezoidal  rule  to  carry  out  the  6 - integrat  ions  in 
Eq.'s  (32),  (33),  and  (34)  and  Romberg  interpolation3,  and  have  used  a 
10-point  Gaussian  quadrature  to  evaluate  the  u-  and  v- integrations. 

For  a value  of  q much  greater  than  two  (2),  there  is  a potential 
accuracy  problem  in  the  evaluation  of  the  integral  in  Eq.  (34).  This  is 
most  easily  studied  if  we  split  the  0- integration  in  Eq.  (34)  in  such  a way 
that  each  region  of  integration  goes  over  only  one-half  cycle  of  the  oscil- 
lating function  cos  (q 0 ) . This  gives  us 


d9  (x,2  t Xj  f - 2Xj  Xj , cos  Qf/6  cos  (q?) 
k n/q 

d9  (x,2  4 Xj  ? - 2 Xj  x,  , cos  0 ) &/b  cos  (q 9 ) . (35) 

(k-  )v/q 

The  factor  (xj2  4 x,  ^ • 2 x(  x, , cos  0)6/s  is  a positive  monotonically- increasing 
function  of  9 in  Eq.  (35),  from  which  we  see  that  the  right-hand  - s ide  of 
Eq.  (35)  consists  of  the  sum  of  a set  of  alternating  sign  terms.  When  there 
is  a large  difference  between  x,  and  xt , , we  note  that  (x,2  + x, ,?  - 2 xt  xs , 

X cos  0)B'°  is  a very  weak  function  of  9 , so  that  the  terms  being  summed 

are  nearly  equal  in  magnitude,  but  of  alternating  sign  --  a situation  that  can 
seriously  stress  the  accuracy  of  the  computed  results. 

To  avoid  this  accuracy  problem,  we  can  expand  the  integrand  in 
Eq.  (34)  in  the  form 

TT  TT 

f d9  (Xj2+  x,  2 - 2 x,  x,  , cos  9F6  cos  (q0)=  ^/6  " d0  (1  - c cos  0)B/e  cos  (q0) 

0 0 

n 

SS/e  d0  £l  - ^ e cos  9 4 ~ j~2~\  ^ cos2  0 + • • • cos  (q0) 

o 

“T5  TT 

= SB/8  n dQ  cos"  0 cos  (q9)  , (36) 

L. . n!  J 

n=  t 0 
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wh  ore 


S = 


+ x, 


(37) 


€ = 2 Xj  Xj , /S 


(38) 


(39) 


Ihe  integral  in  the  final  form  of  the  right-hand-side  of  Eq.  (36)  can  be 
shown  to  have  the  value4 


( 


d8  cos"  9 cos  (qQ)  =J  o 


if  n < q 
if  n-t-q  = odd 


. (40) 


r n ! 


^Ti)'f7L)i 


if  n-tq  = even 


so  that  Eq.  (36)  may  be  rewritten 


a s 


dS  (x,  + x,  _ -2xtXj,  cos  Q)5/  cos  (qQ) 

(e/2)*+** 


= 77  gR/6 


l , 

, -o  (q+ni ) ! m ! 


(41) 


For  «<.5  . which  corresponds  ..  x,  and  x,.  signif  ic.ntly  d ifferenl 

air,..  Eq.  (41)  is  fairly  rapidly  convergent  and  a 40-term  summation 

yields  sufficient  accuracy.  For  « > . 5 . x,  and  x,.  ere  .nfflcl.ntl, 

close  in  value  that  (x  3 -t  x 3-2x  v rno  ms/«  • 

• t'  ^ xi  xi'  cos  ls  a significant  function 

of  9 . and  the  evaluation  of  the  integral  in  Eq.  (34)  can  proceed  by 

straightforward  numerical  quadrature  without  excessive  loss  of  accuracy. 
For  q greater  than  two  (2),  we  have  evaluated  the  0 -integration  in  Eq.  (3 
using  either  ordinary  nun  erical  integration  techniques  if  2 x,  x,  , /(x3  4 Xj  / 
is  greater  than  one-half,  or  used  Eq.  (41)  for  values  less  than  one-half! 
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With  the st  numerical  techniques,  the  Kqs  (i,  i ')  matrices  were 
evaluated  lor  q 0 to  41  . These  results  are  listed  in  Table  I,  giving 

he  matrix  in  upper  right  triangular  form.  With  these  matrices,  we  were 
then  able  to  carry  out  a determination  of  the  eigenfunctions  and  eigen- 
values. These  a*-e  listed  in  Table  II.  The  eigenfunctions  as  listed  here 
have  been  restore!  to  their  unsymmetrized  form,  i.  e.  , UpMx)  instead 
of  (x)  by  making  use  of  Eq.  (26). 

The  set  of  all  eigenvalues  were  rank-ordered  without  regard  to 
p-  and  q -values.  In  this  procedure,  we  counted  each  eigenvalue  twice 
if  its  q -value  was  not  zero,  since  it  then  applied  to  both  q and  -q 
This  set  of  eigenvalues  is  listed  in  Table  III.  We  note  that  the  leading 
eigenvalue  for  each  value  of  q for  q ^ 4 appear  in  order  according  to 
the  \alue  of  q . Since  we  only  worked  with  q ^ 41  , and  since  the  lead- 
ing e i genva  lue  for  q 41  is  the  569  , we  can  probably  consider  the  list 

of  eigenvalues  complete  up  to  the  569,  or  thereabout.  The  sum  of  all  the 
eigenvalues  listed,  of  which  there  are  1660,  is  0,10512*  . Th is  is  in 

good  agreement  with  the  value  of  0.  1056  expected  for  the  total  of  all 
eigenvalues,  as  derived  from  an  earlier  work  which  considered  the  ex- 
pected mean  square  wavefront  distortion,"  We  note  that  the  cumulative 
sum  at  the  569tl  eigenvalue  is  0.104708/0.  10527  = 99.60%  of  the  total  of 
the  eigenvalues  listed,  and  0.  104708/.  1056  = 99.  15%  of  the  total  of  all 
the  eigenvalues. 

With  this  list  of  eigenvalues,  it  is  possible  to  proceed  immediately 
to  the  imaging  probability  evaluation  aspect  of  the  problem.  This  we  take 
up  in  the  next  section.  Before  turning  to  that,  however,  we  first  note  that 
because  of  the  ease  with  which  we  < ould  adapt  our  mathematics  to  the  case 
in  which  tilt  is  considered  to  be  a significant  wavefront  distortion,  we  have 
carried  out  such  calculations.  As  noted  previously,  this  involves  nothing 
more  complex  than  using  Eq.  (34)  in  place  of  Eq.  (33)  for  q 1 . With 


- 38  - 


this  replacement,  we  obtain  the  symmetrized  kernel  shown  in  Table  la, 
and  the  eigenvalues  and  eigenfunctions  {with  symmetri zation  removed) 
shown  in  Table  Ha.  It  is  particularly  interesting  to  not-  that  the  first 
eigenvalue  in  this  case  is  about  3g  times  larger  than  the  sum  of  all  the 
eigenvalues  when  tilt  effects  are  suppressed  as  not  contributing  to  effec- 
tive wavefront  distortion,  and  that  the  first  eigenfunction  appears  to  be 
very  nearly  a simple  tilt.  We  also  note  that  the  second  and  subsequent 
eigenvalues  with  tilt  distortion  allowed  (Table  Ila)  are  very  nearly  equiv- 
alent to  the  first  and  subsequent  eigenvalues  when  tilt  distortion  is  not 
allowed  (Table  II,  q 1 ). 

Probability  Integrals 

Having  the  list  of  eigenvalues  given  in  Table  IH  reliable  out  to  the 
t h 

569  eigenvalue,  and  thus  reliably  containing  more  than  99%  of  the  sum 
of  all  eigenvalues,  we  are  now  in  a position  to  start  the  evaluation  of  the 
probability  integral  governing  the  performance  of  a CENSORING  system. 
We  recall  that  this  integral  is  given  by  Eq.  (20).  Rewriting  this  to  work 

with  the  n-notation  (overall  rank  order)  of  Table  III,  rather  than  with  the 
p,  q-notation,  we  write 


CENSOR 


(2  n c„sr1/a 


! 

Sp  h 


dxn  exp  (-£  x,2/c  2) 


(42) 


where  the  "Sph  limit  on  the  n -d ime ns iona  1 integration  corresponds  to 
the  constra  int 


) 

B 


(43) 


In  ai  cordance  with  Eq.  (19),  with  denoting  the  n1*1  eigenvalue  in 
Table  III,  we  h ave 


a 


2 

n 


(44) 
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for  an  aperture  of  diameter  D with  wavefront  distortion  characterized 
by  the  length  r0 

We  plan  to  evaluate  the  n-dimens ional  integral  in  Eq.  (42)  by 
Monte  Carlo  techniques.  As  an  immediate  reduction  in  the  magnitude 
of  the  problem,  we  recognize  that  0B2  decreases  rather  rapidly  with 
increasing  n , and  that  beyond  some  cut-off  value  of  n , the  expected 
value  of  the  sum  of  the  square  of  the  random  variables  is  very  tightly 
constrained  to  the  sum  of  the  , with  very  little  variability.  The  real 

variability  in  the  overall  sum  of  the  squares  comes  from  the  much  fewer 
x,'s  for  which  jn2  is  large  and  decreases  rapidly  with  increasing  n . 

To  establish  the  cut-off  value  of  n namely,  N , we  arbitrar- 
ily allow  the  random  variables  xB  beyond  the  cut-off  to  have  an  expected 
value  of  0.  1 for  the  sum  of  their  values  squared,  i.  e.  , 

00  oc 

< X,2)  = a,2  = °.  1 . (45) 

. — 

\ Nc 

Since  we  know  that 

X 

) a 2 = 0.  1345  (D/r0JP/s  , (46) 

L " 

i 

then  we  can  obtain  N from  the  running  cumulative  value  of  the  eigen- 
values in  Table  III.  We  write  in  accordance  with  Eq.  (44) 

Ne 

-(D/r„)Ey3  * 2 = 0.  1345  (D/n.jf5'3  - 0.  1 , (47) 

TT  ° _ n 

1 

from  which  it  follows  that 

\ 

) « 2 = 0.  1056  - 0.  078o4/(D/ro)^/3  . (48) 

L * 
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At*  see  that  to  avoid  pushing  the  cut-off,  N , beyond  the  reliable  limit 
ol  our  eigemalue  table,  i.  e.  , beyond  about  569,  the  largest  value  of  D/r0 
we  tan  use  is  14.60  (which  we  take  to  be  15). 

For  each  value  of  D/ r0  , we  determine  N.  in  accordance  with 
Fq.  (48)  and  Table  II,  and  then  proceed  with  the  evaluation  of  the  truncated 
probability  integral  of  Eq.  's  (42)  and  (43),  which  we  now  rewrite  as 

P-'nso,  = l (2"  cnr)"1,‘  ’ d.^  exp  (-£  \z/cn2)  , (49) 

where  now  the  limit  on  the  N -dimensional  integration  is  given  by  the 
const  raint 


In  the  evaluation  of  the  integral  in  Eq.  (49)  by  Monte  Carlo  methods , 
there  are  a number  of  approaches  to  the  random  sampling  that  can  be 
used.  First,  and  most  obvious,  we  consider  selecting  points  uniformly 
distributed  in  the  hypersphere  of  Eq.  (50),  and  evaluate  the  integrand  of 
each  point.  Unfortunately,  the  integrand  will  be  very  small  for  most  of 
the  points  selected,  since  many  values  of  ~n  , for  large  n , will  be  much 
less  than  unity.  This  will  give  very  poor  sampling  efficiency  and  an  unman- 
ageably large  number  of  samples  will  be  needed  to  yield  even  modest  accur- 
acy, 

A second  approach  to  the  random  sampling  is  to  select  the  random 
points  in  accordance  with  the  gaussian  probability  distributions  for  each 
dimension  inherent  in  the  integrand  in  Eq.  (49).  Then  the  integral  would 
be  evaluated  by  counting  a one  for  each  such  randomly  selected  hyper-space 
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point  that  satisfied  Eq.  (50),  and  zero  for  each  point  that  did  not,  and 
taking  the  average  of  these  counts  of  one  and  zero.  Unfortunately,  this 
method  also  suffers  from  the  problem  that  an  unmanageably  large  number 
of  sample  points  are  needed  to  obtain  acceptable  accuracy  in  the  integral 
evaluation.  In  this  case,  the  problem  is  associated  with  the  smaller  values 
of  n , especially  for  larger  values  of  D/r0  . The  values  of  a.  are  so 
large  that  the  gaussian  distribution  of  x„  causes  most  points  selected  to 
lie  outside  the  hypersphere  defined  by  Eq.  (50). 

To  get  around  the  problems  of  both  the  first  and  the  second  methods 
of  sampling,  we  used  an  arbitrarily  chosen  sampling  distribution  Tl,  (xB) 
for  each  of  the  Nc  dimensions  of  the  integral.  To  compensate  for  this 
method  of  choosing  the  samples,  it  is  merely  necessary  to  introduce  a 
factor  of  1 '/'j  ^(xJJ  into  the  integrand.  We  choose  (xt)  to  match 
the  gaussian  distribution  with  variance  for  the  larger  values  of  n in 

the  integration,  for  which  ca:  is  small  and  there  is,  therefore,  no  ten- 

dency to  pick  values  of  x,,  that  are  incompatible  with  Eq.  (50).  For  the 
smaller  values  of  n , with  larger  corresponding  values  of  , we 

chose  a gaussian  distribution  with  a variance  302  . Here  eg  is  the  same 

for  all  of  the  dimensions,  and  significantly  less  than  the  corresponding  eg 
values.  To  establish  the  transition  between  what  we  have  called  the  large 
values  of  n and  the  small  values  of  n , we  determined  a transition  value, 
Nr  , which  would  satisfy  the  requirement  that 


where 
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0.  9 


(51) 


(52) 


The  value  of  is  directly  obtainable  from  the  eigenvalues  of  Table  III, 

using  Eq.  (44). 
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Using  a gaussian  sampling  distribution  with  variance  oQz  = a^2 
for  variables  x,  , Xj  , Xj  , . . . xN  , and  variance  an2  for  the  variables 
x*,  + j . \ +a  - • • • - xu  > vve  have  found  that  Eq.  (50)  is  satisfied  between 
one -third  and  two-thirds  of  the  time  by  the  hyper -space  random  vector  so 
chosen.  Using  this  sampling  procedure,  and  noting  that  now  the  probability 
integral  of  Eq.  (49)  has  the  form 
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The  evaluation  of  PCENS0R  was  carried  out  in  accordance  with  Eq.  (53) 
using  samples  of  100  points  for  various  values  of  D/ rQ  . By  repeating  the 
evaluation  a number  of  times,  it  was  possible  to  obtain  an  estimate  not  only 
of  P,ENs:,R  , but  also  of  the  rms  uncertainty  in  our  answer.  In  Table  IV, 
we  list  our  results.  These  results  are  the  basic  objective  of  our  numerical 
exercise.  In  the  next  section,  we  discuss  the  interpretation  of  these  results. 


Discussion  of  Results 

The  basic  results  are  those  presented  in  Table  IV  and  we  shall  center 
our  discussion  about  these  results.  The  very  large  variation  of  PCEN30R  with 
D/r0  is  apparent  from  even  a cursory  examination.  In  Fig.  1,  we  have 
plotted  PrEN30R  as  a function  of  ( D / rQ )2  . It  is  interesting  to  note  how  well 

the  data  is  fit  by  an  exponential  dependence  on  aperture  area.  This  is  in 
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good  agreement  with  an  earlier  conjecture  by  Hufnagel,3  though  the  coef- 
ficients of  the  fit  are  significantly  different  from  those  suggested  by 
Hufnagel.  We  find  that  the  data  is  well  represented  by  the  relationship 

Pc£ns0r=  5.  6 exp  [-  0.  1 557  (D/rr)*]  . (55) 

at  least  for  values  of  D/ rQ  s 5 . 

It  is  interesting  to  remark  that  if  a CENSORING  experiment  were 
performed  with  D/rQ  = 1 5 , as  would  be  the  nominal  condition  for  a 1 . 5 m 
telescope  (with  r0  nominally  equal  to  0.  1 m ),  then  the  probability  of 
getting  a good  picture  in  a single  short  exposure  would  be  about  3.  4 X 10“  15  . 
If  independently  distorted  wavefront  short  exposures  could  be  obtained  at  the 
rate  of  100  per  second,  it  would  take  more  than  800  million  hours  "on  an 
average"  to  get  a good  picture,  i.  e.,  one  for  which  the  average  wavefront 
distortion  over  the  aperture  was  less  than  one-radian.  If  the  aperture  dia- 
meter were  reduced  to  1 m , so  that  D/r0  = 10  , the  probability  would 
be  about  l.  1 X 1O-0  , and  the  expected  waiting  time  to  get  a good  picture 

would  become  about  2.  5 hours  (if  we  can  get  100  indepenuent  wavefront 
distortion  samples  per  second).  With  a 0.  7 m diameter  aperture,  the 
waiting  time  shrinks  to  only  3.  5 seconds.  Clearly,  in  a CENSORING  ex- 
periment, it  is  critical  to  know  what  rQ  is  and  to  not  make  the  aperture 
diameter  much  larger  than  about  7 r0  , unless  very  long  waiting  times 
are  acceptable. 

We  note  that  in  certain  cases,  astronomical  seeing  with  rQ  values 
in  excess  of  0.  1 5 m have  been  reported.7  In  such  cases,  it  would  be  quite 
appropriate  to  attempt  a CENSORING  experiment  with  aim  diameter  aper- 
ture, but  it  is  critical  that  the  aperture  be  properly  stopped,  and  this  re- 
quires current  knowledge  of  rQ  , and  appropriate  planning  in  the  imple- 
mentation of  the  experiment. 


44 


C tNSCR 


I 


icr'2 


10-1* 


Figure  1.  Probability  of  Obtaining  a Good  Short  Exposure  Image  as  a 
Function  of  Aperture  Diameter.  A good  image  is  defined  as 
one  with  less  than  one  radian2  effective  wavefront  error 
(i.  e.  , wavefront  error  excluding  tilt)  over  the  aperture. 
Aperture  diameter  D Is  measured  in  units  of  the  wavefront 
distortion  length,  rQ  . PttN90l,  is  the  probability. 


Table  III 
Eigenvalue  List 


Each  eigenvalue  and  its  associated  values  of 
p and  q are  listed,  for  all  values  covered 
in  Table  II.  When  q t 0 , the  eigenvalue  is 
considered  to  be  listed  twice,  as  indicated  by 
the  nature  of  the  overall  rank-order  column, 
N , and  by  the  Q ,alue  column. 
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13 

11,-11 

9 

20,-20 

14 

J.  -9 

6 

29.-20 

5 

31  ,-31 

16 

5,  -5 

15 

7,  -7 

3 

39,-38 

10 

19,-18 

9 

23,-23 

17 

4,  -4 

11 

16,-16 

12 

14,-14 

4 

35,-35 

1 3 

12.-12 

14 

10,-10 

7 

26,-26 

15 

8,  -8 

IB 

3,  -3 

19 

2.  -2 

16 

6,  -6 

9 

21.-21 

5 

32,-32 

6 

29,-29 

3 

39,-39 

10 

19,-19 
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78? 

783 

.000000  719 

784 

785 

.000000737 

78ft 

787 

.000000715 

788 

789 

.000000715 

790 

791 

.000000733 

792 

793 

.000000731 

79* 

795 

.000000728 

79ft 

797 

.000000728 

798 

799 

.000000728 

800 

801 

.000000720 

flo? 

8o  3 

.000000716 

HO* 

805 

.000000715 

80. 

807 

.000000712 

808 

809 

.000000706 

810 

81  1 

.000000706 

81? 

813 

.000000705 

81* 

815 

.000000  703 

HI  ft 

81  7 

.000000597 

818 

81  9 

.00000069 7 

820 

8?1 

. 000000691 

82? 

H?3 

.000000689 

82* 

H?6 

.000000688 

B?6 

8?  7 

.000000687 

8?R 

H?9 

.00000068  7 

830 

811 

.000000673 

«12 

813 

.000000670 

81* 

815 

.000000668 

«3ft 

817 

• 0000006ft  1 

818 

819 

. OQ0000661 

8*0 

8*1 

.000000656 

8*2 

8*3 

.000000656 

8*  * 

8*5 

,0000006*9 

8*  ft 

8*  7 

. O00000‘>*6 

8*8 

8*9 

.0000006*5 

850 

«S1 

. 0000006*5 

852 

853 

. 0000006** 

85* 

855 

.0000006** 

85ft 

85  7 

. 000000^*3 

858 

859 

.0000006*2 

BhO 

8ft  1 

.0000006*0 

8ft? 

Hft3 

.OO00O0611 

86* 

8ft5 

,000000611 

86ft 

8 ft  7 

.000000628 

BbH 

869 

.000000628 

870 

871 

.000000620 

87? 

873 

.00000061 9 

h 7* 

871 

.000000617 

87ft 

87  7 

.000000610 

878 

8 7 9 

.000000606 

8b  0 

8 H 1 

.000000605 

r.u*i,  sjm.  p o 


.io*9i  no? 

1 r 

6. 

-5 

. 1 0*91  2777 

16 

7, 

-7 

. 1 0*91 *?*7 

15 

9. 

-9 

.10*915717 

8 

24, 

.24 

.10*91 71 8* 

11 

17, 

-17 

. 10*91 86*6 

14 

lit 

-11 

. 1 0*9201 03 

12 

Iftt 

• 15 

.104121559 

13 

13, 

-13 

.10*923015 

18 

4, 

-4 

. 10*92**55 

4 

36, 

-36 

. 1 0*925886 

7 

27, 

-27 

.10*927116 

17 

6 1 

-6 

. 1 0*9?87*0 

9 

22, 

-22 

. 10*9301  53 

6 

30, 

-30 

. 1 0*911565 

5 

33, 

-33 

. 1 0*912975 

16 

-8 

. 10*91*181 

3 

40, 

-40 

.10*935775 

10 

20, 

-20 

. 10*9  17] 69 

15 

10, 

-10 

. 10*93855? 

14 

12, 

-12 

, 1 0*939911 

1 1 

1 8, 

-1  8 

. 10*9*1106 

1 3 

1*, 

-1* 

. 1 0*9*2ft80 

8 

26, 

-25 

. 1 0*9**053 

12 

1 6 , 

-16 

. 10*9*5199 

* 

37, 

-37 

. 1 0*946719 

7 

28, 

-28 

. 10*9*807* 

9 

21, 

-23 

, 1 0*9*9197 

5 

3*, 

-34 

. 1 0*950718 

ft 

31  , 

-31 

. 1 0*952030 

3 

41  , 

-41 

. 1 0*9531*2 

10 

21  , 

-21 

. 10*95*619 

1 l 

19, 

-19 

. 1 0*955911 

\ 2 

17, 

-17 

, 1 0*057220 

8 

26t 

-26 

. 10*958509 

1 3 

15, 

-15 

, 10*959797 

14 

13, 

-13 

, 1 0*96108* 

1 5 

n , 

-1  1 

. 1 0*962171 

16 

9, 

-9 

. 10*961655 

1 r 

7, 

-7 

. 1 0*9ft*9  35 

1« 

5, 

-5 

. 1 0 *Oftft 1 97 

4 

38, 

.38 

. 1 0*96  7*59 

1 9 

3, 

-3 

. 1 0*968715 

7 

29, 

-29 

. 1 0*9ft9970 

9 

2 *,-?* 

. 1 0*971210 

ft 

12,-32 

, 1 0*972**8 

5 

35,. 

-35 

. 1 0*97 3 ftH? 

1 o 

22, 

•22 

. 1 0*97*901 

1 1 

2o , 

-20 

. 1049  761  1 * 

H 

27,-27 

.10*977123 

1? 

1 8 , - 1 8 
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riiM.  s.im. 

P 

9 

88?. 80  H 

. noonon-iQQ 

. 1 049785? 1 

1 3 

1 5 , -1  b 

884. 80S 

.OOOUOO-.92 

. 1 04979705 

4 

39,-39 

H0b*8R  / 

. 3 00  no 0 591 

. 1 04980087 

14 

14,-14 

808,889 

. 0 00(1 00  HQ  1 

. 1 04QO?0b9 

9 

25,-?? 

890 . 89 1 

. Oil  000059) 

. 1 0498  3?5 1 

7 

30.-30 

89? 

. noooooiOB 

. 1 0490  3839 

?n 

0 

091,894 

, on  00 MO  HR  3 

. 1 ('49850  04 

5 

3b. -3h 

895.89b 

.00000070? 

, 1 0 4 9 Ob  1 59 

b 

33,-33 

89  7 .898 

.000000501 

. 1 0498  7331 

IS 

1?,-1? 

HQ  -J  « y 0 u 

. 000000-.H  1 

. 1 0498H49? 

1 0 

23. -?3 

Q li  1 .902 

. ) 1*0  0 005  7? 

. 1 04Q89blh 

1 1 

?1 .-?) 

9li  3 * Qn** 

, 000(1  005  71 

. 1 04991(778 

H 

2h. -?H 

Q(i  ->  » Q Ob 

. 000000  -.bb 

.10409191 1 

1b 

10,-K) 

90  7.908 

. 0 00ii00-»b3 

. 1 04993037 

1? 

19,-19 

9O9.910 

,000000557 

. 1 040941  'I 

4 

40,-40 

9 1 1 » 9 1 2 

. OOOlMlO-ib  7 

, 1 0 4 99b? 5- 

9 

?b,-2b 

91  3.91  <♦ 

. OOOOOOibh 

. 1 04995378 

7 

31  .-31 

915.91b 

. OOOOUObS? 

, 1 04QQ  74O 1 

13 

1 7 » - 1 7 

91 7.91 8 

. 000000--4Q 

, 1 04998070 

b 

34.-34 

91 9,9? 0 

. 00000054H 

, 1 04999573 

5 

3 7,-3  7 

9? 1 ,9?2 

. O00d0054b 

. 1 ObOOOTbb 

1 0 

?4,-?4 

923. 9?4 

. 00000054  3 

, 1 ObO  0 1 05 1 

1 7 

8,  -8 

9 ? 5 « 9 ? h 

. 0 000  00 ->4  1 

. 1 050  0 ?9  3? 

?() 

1 . “1 

9 ? 7 , Q ? H 

. 00000053H 

. 1 05()04008 

8 

29. -?9 

9 ? 9 , 9 1 0 

.onoo 00537 

, 1 0500500? 

1 4 

15,-15 

931  .932 

. 000000535 

. 1 0500515? 

1 1 

2?.-22 

933,934 

.0000005?? 

. 1 (1500  7?03 

7 

3?. -3? 

9 3 5.93b 

. 0 0 0 0 0 0 b ? 5 

, 1 Ob008?5? 

9 

27,-27 

93  7 .939 

« 0 0 0 0 (.  0 5 ? b 

. 1 05009301 

4 

41,-41 

939.940 

.0000005?? 

,105010344 

1? 

20.-20 

94  1 » -*42 

.00000051 7 
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is 
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94  3,944 
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,1 0501 2414 

6 

35,-35 

945 » 94b 

.('0000051  7 
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S 

38,-38 
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..30000051  ? 
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10 

?5 , -25 

949.9mi 

,00000050b 

.105015489 

8 

30,-30 

9b) . 952 

.000000 50 b 

. 1 ObO 1 b500 

1 >3 

1 0 , - 1 B 

95  3.954 

. 000000^01 

,10501  750? 

1 b 

b , — b 

955,95b 

,000000498 

, 1 050  1 8498 

1 1 

23,-23 

957,958 

, Ou()0004Qb 

, 1 050 1 9490 

7 

33.-33 

959.9b(J 

,000000494 

, 1 OSO?0479 

9 

?o , -28 

9b 1 » 9b? 

. 3000004H9 

. 10n0?l457 

lb 

n ,-n 

9 b 3 , 9 b 4 

,000000459 

, ) 050?24  3b 

b 
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,000000488 

, 1 050?  3411 

5 

39,-39 

9b  7 , Qbrt 

, 000000454 

, 1 05024380 

14 

lb, -lb 

9b9. 370 

. O0O0004H? 

, 1 05025343 

1 ? 

21 ,-?l 

971  .97? 

,000000450 

. 1 050?b303 

1 0 

?b,-2b 

973,974 

,((00000490 

. 1 050?  7?b? 

8 

31 ,-31 

97b»97b 

, " Q 0 0 0 0459 

, 1 050?8?01 

7 

34,-34 

977. 97H 

,00000045b 

, 1 f(b0?91  3? 

9 

29,-29 

979.980 

. 000000453 

, ) 05030050 

1 1 

?4 , - ?4 
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9H1 

)«? 

963 

■in* 

96-, 

966 

96  7 

9p  j 

39  0 

99  1 

wq  ? 

99  3 

994 

995 

996 

99  7 

99H 

99  9 

1 nun 

1001 

1(H)? 

10  0 3 

100  4 

1005 

lOO*' 

1007 

1 Ol)M 

1 009 

lOlii 

mu 

1 0 1 ? 

1013 

1 0 1 4 

10  14 

1 0 1 6 

1017 

1 Ol  4 

1019 

1 p?0 

1021 

1 o2? 

10?  1 

1 (t2<* 

1 U?4 

1 0?6 

10?  7 

1 024 

1 0?  9 

loin 

1031 

1 nJ? 

1033 

1 o3<* 

1034 

1 P 34 

1037 

103m 

1 03  9 

1 p 4 o 

1041 

l P4? 

1 04  3 

In44 

104  5 

1 P4<1 

1 047 

1 046 

1049 

1 n5o 

1051 

1 o5? 

1053 

1 P54 

1 054 

1 056 

1057 

1 pSh 

1 049 

1 061) 

1061 

1 06? 

1063 

1 0 6 4 

1064 

1 064 

1 06  7 

1 068 

1 069 

1 o 70 

1071 

107? 

1073 

10  74 

1074 

1 0 74 

107  7 

1 0 76 

1 079 

1 OHO 

.O0f)000<*63 
.fioonoo*.^? 
. nooooruM 
.000000<.S7 
. nooooo<»^3 

. <H/O0  00<*4  9 
.01)0000447 
.000000^44 
. 000000444 

• 0OOt,00439 
. o o n o o o 4 3 7 
.nnnonn'.  17 
. O(»nod04  14 
.000000430 
. 000u004?fl 

• o n n o o o -*  ? i 
, n(.000  T4P0 
, o ft  ft  ftt)  o 4 1 9 
. 0 0 0 n 0 0 4 1 9 

.noouonoh 

. 0(100004  1 1 
.000000410 

.oooiiuo4  0 8 

. OO00004O4 
.000000401 
.000000399 
.000000398 
. o 0 0 0 0 0 3 9 S 
. 00000039? 
.O0O0003Q0 
.000009  1«6 
, O0O00O3H1 
.00^000301 
. 000000  17H 
. 000000373 
. OGOOOO  175 
. '>00000375 
.000000370 
. 000000344 
. O0000O  145 

.O0O000343 
. 000000  140 
. 000000  158 
.000000357 
. 900000  15? 
.000000349 
.000000346 
. 000000344 
. 000000  * 4 J 
. OUOoOO  14? 


ft'*.  Sum. 

. 1 05030944 
.105031904 
. 1 05032430 
. 1 05033743 
. 1 05034450 
. 1 05035547 
. 1 05036442 
.105037330 
.105038217 
. 10503^094 
. 105039969 
. 1 05040443 
.105041713 
. 105042573 
. 1 05043429 
. 1 05044270 
. 1 050451 1 0 
. 1 05046949 
. 105046746 
. 1 05047M  7 
.105048439 
. 1 05049260 
. 1 05050074 
, 1 05050444 
.1 05051444 
.105052445 
, 1 05053241 
.105054071 
. 1 OS05445S 
.105055434 
, 1 05056406 
.105057169 
.105057931 
.105058644 
,105059443 
. 1 050401 93 
. 1 05060943 
. 1 05061642 
, 1 05062414 
, 1 050631 44 
.105063469 
, 1 05064549 
. 1 05065304 
.105046014 
.1050667?? 
. 1 0506  7420 
.1050681 12 
. 1 0504«»01 
. 1 05069446 
. 1 050701  69 


P 0 


6 

37, 

-37 

5 

40, 

-40 

13 

19, 

-19 

15 

u, 

-14 

4 

32, 

-32 

10 

27, 

-27 

17 

9, 

-9 

7 

34, 

-35 

12 

22, 

-22 

6 

38, 

-38 

9 

30, 

-30 

5 

41  , 

-41 

14 

17, 

-17 

1 1 

25, 

-25 

8 

33. 

-33 

7 

34, 

-36 

1 0 

28, 

-28 

13 

2o, 

-?0 

15 

12. 

-12 

6 

39, 

-39 

9 

31  , 

-31 

19 

4, 

-4 

12 

23, 

-23 

8 

34, 

-34 

15 

15. 

-15 

1 1 

26. 

-26 

7 

37, 

-37 

6 

40. 

-40 

10 

?9, 

-29 

14 

18, 

-18 

9 

32. 

-3? 

8 

35, 

-35 

13 

21. 

-21 

18 

7, 

-7 

7 

38, 

-38 

6 

41  , 

-41 

12 

?4, 

-24 

11 

27, 

-27 

10 

30, 

-30 

1 7 

lo, 

-10 

9 

33. 

-33 

8 

36, 

-36 

7 

39, 

-39 

16 

13. 

-13 

15 

16, 

-16 

m 

19, 

-19 

13 

22. 

-22 

1? 

25, 

-25 

1 1 

PP. 

-28 

1 0 

31  , 

- 3 1 
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£1(>EN 


ru*.  sum 


0 


1 


1081  . 

1080 

1083. 

1 084 

1085. 

1 0M8 

1087. 

1 08H 

1089* 

1090 

1091  . 

1090 

1091. 

1 09fc 

1 09s  • 

1 098 

109  '» 

1 098 

1.(99. 

1 100 

1101  . 

1 1 00 

1103. 

1 104 

1105. 

1 1 08 

110  7. 

1 1 OH 

1109. 

1 1 10 

1111. 

1 1 10 

mi* 

HU 

1115. 

1 1 1 8 

1117. 

1 1 18 

1119. 

1 12L 

1101* 

1120 

11  ? 3 , 

1 1 04 

110*. 

1 1 28 

1107, 

1120 

1 109. 

1 1 JO 

1131. 

1 1 30 

1131. 

1 1 34 

1135, 

1 1 J8 

113/. 

1 1 JH 

1139. 

114" 

1141. 

11*0 

1U3. 

1 1 44 

1 1 49. 

1148 

1147. 

1 1 

1 l 4 

l 1 SO 

1 151. 

1 1 50 

11M. 

1 1 5<» 

1 1 58, 

1 1 58 

1157. 

1 1 58 

1189, 

1 1 80 

118)  . 

1 1 80 

1183. 

1 184 

1148. 

1 1 8» 

1187. 

1 1 8t 

1189. 

1 1 70 

11/1. 

1 1 7? 

117  1. 

H/4 

1 1 78. 

1 1 7 4 

1177. 

1 1 7 u 

1 1 79. 

I 1 8<i 

.000000741 
. 000000140 
. 000000339 
. 000000301 
. 00OO0O371 

.oooooo i?r 

.ooOOOO  119 
.000000 118 
.000000 118 
.000000115 

. 0 U o 0 0 o 3 1 ? 

.OOOOOO  109 
. OOOuOO 
.000000103 
. OOOOOO 10  1 
. (1 0 0 0 O 0 0 9 H 

.000000097 

.00000 0 095 
. 0 (I  o 0 0 0 0 9 1 
.000000068 
. OOOOQO0R8 
.000000095 
.000  0 00 ->81 
.000000090 

.0001100079 
,O00O0O?74 
. ooouoop r? 
. 000000->7o 
. O0O(>00?88 
. no  Onoo  088 
. oooooo-' *>n 
."0000008(1 
.000000880 
. 0 0 0 r.  o o ? 5 8 
. 0(10000854 
,Oi,0o0O851 
. 00P00O?80 
. 0001)0084  7 
. 0 0 O K)  0 8 4 <* 
.0(1  0 00  0 84  . 
.000000899 
. 0 (1 0 0 0 0 8 1 8 
.00 00 00838 
.OOOf  00  8 98 
. OH0000888 
.OO(K»OO>?0 
. oooo jn??b 

. ')00(i0fl8?4 
, 00000080? 
. OOOQOO00O 


. 1 05070081 
. 1 05071530 
. 1 05072009 
. 1 0507095? 
. 1 05073404 
.105074133 
. 1 05074770 
. 1 08075407 
. 1 05078040 
. 1 05078889 
. 1 050  7 7094 
.105077910 
.105078501 
. 1 050791 08 
. 1 05079798 
. 1 050«03?<» 
. 1 0500091 8 
. 1 08001500 
. ) 05000090 
. 1 05000883 
. 1 05003035 
. 1 0508  1905 
. 1 08094370 
. 1 0S0B4O30 
. 1 05005487 
. 1 0 80080  35 
. 1 05008570 
.105007119 
.10809/854 
. 1 05000107 
. ) (15000708 
. 1 08009009 
.105009749 
. 1 05090085 
. 1 05090774 
.105091078 
. 1 05091 777 

. 1 08090070 

. 1 PSPO075W 
. 1 05093045 
. 1 0509  1704 
. 1 05094000 
. 1 0509487? 
.1 080951  44 
, 1 05095*0  1 
. 1 05098057 
.105098509 
. 1 05098958 
. 1 08(197400 
. 1 0809/041 


9 

34 

.-34 

R 

37 

.-37 

7 

40 

,-40 

7 

41 

.-41 

8 

J« 

. -38 

9 

35 

.-35 

10 

3? 
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29 

.-29 

12 

26 

13 

23 
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14 

2o 
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15 

17 
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16 

14 
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8 

39 
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9 

38 

,-38 

1 0 

33 
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1 1 

30 

,-30 
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.-07 

R 

40 

,-4  0 

13 

04 

.-24 

1 8 

9 

* • H 

9 

37 

.-37 

14 

01 

.-21 

1 0 

34 

.-34 

1 1 

31 

. -3 1 

18 

1« 

• - 1 W 
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41 

.-41 

12 

29 

, — 0R 

9 

35 

.-38 

13 

05 

.-05 

1 0 

14 

.-38 

16 

16 

.-15 

1 9 

8 

. -5 

1 1 

30 

.-30 

14 

0? 

.-?? 

9 

)« 

• — 99 

10 

09 

.-09 

1 0 

34 

, —38 

1 7 

1 ? 

.-1  ? 

18 

19 

.-19 

1 3 

?4 

. —06 

1 1 

33 

.-3  3 

9 

40 

. -*  0 

1 0 

37 

.-3  7 

12 

30 

.-10 

14 

23 

.-03 

16 

1 4 

,-16 

9 

41 

.-41 

1 1 

34 

• - .14 

. 58  . 


9 


1161.1182 
1 1 6?»  ) 
i i*5.  i i of- 
Ilf-  1%  1 1 88 
I 169.  1 1 90 
1191*11 

) 105.1  1 9#. 
119  7*11  9** 
1199*1 ? 0 0 
l if  0 1 .1202 
1203.120* 
1 2<)5«  1 20* 
1 2c  7. 1 ?0» 
12), 9.1210 
1 <?  1 1 . 121* 
121  M?U 
12)3.1216 
121  r. i p i ft 
121  9. 1 ??0 
1221.1 22* 
1223.1 22* 
122*. 122* 

1 2?  7 , 1 ?2* 
1229.  |?J(i 
1231 .1232 

1239.1 234 

1235.1235 
1217.1 23* 

1 239.  | ?*n 
12*1*1 2*2 
12*1.1 2** 

1 2*5. 1 2** 
12*  7. 1 2** 
12*9.1 2b" 
1251.1 272 
1251.12b* 
1255.1225 
1 c.5  7.  ) 2b6 

l?59»)2bO 

12M.126? 

1253.126* 

1 255.1 255 
1257.1265 
1259.1270 
1271.J272 
1273.127* 
1275. 1275 
12/7. 127* 
1279. 1280 


VALUE 

."0O.;0O21  7 
.00000021 7 
.00000021* 
,000000212 
.00000021 o 
. 000000*05 
.000000*0* 
.000000201 
.000000*01 
.000000199 
.000000195 
. nyOilOu  1 93 
.0000UO1 92 

• ooonoo i a 9 
.OOOOOO 1 5« 
. O0OQOO 1 5* 
• O0O000 1 °? 
.OOO0001 50 
.000O001 79 
.000000)75 
.000000157 
.000000 157 
.0U0J00157 
.000000)57 
.000000157 
.000000)57 
.OQO0CO157 
.000000157 
.000000157 
.000000156 
.000000)57 
.000000155 
. O00O00 1 5* 
.000000)52 
.0000001*9 
.0000001*7 
• OOOuOOl *5 
.0000001*5 
.0000001*2 
.0000001 *0 
.000000139 
. "0O00O 1 35 
.0000001 35 
.00000013* 
.000000131 
.00000013) 
.000000)25 
.000000127 
.000000)26 
.000000)2 2 


ru*.  Sin*. 

, 1 050O627* 
. 1 05005705 
. 1 050991 35 
. 1 05009550 
. 1 05099950 
. 1 051 00391 
.105100799 
. 1 051 *1201 
.105101503 
. 1 051 02000 
. 1 05102369 
.105102775 
. 1 051 031 55 
.105103536 
.105103912 

. 1 0510*26) 

. 1 051 0*5*5 

.105)05005 

. 1 05105353 

.105105716 

.105105053 

. 1 05J  05366 

.105106722 

. 1 05107057 

•1051* 7392 

.105107725 

. 1 05108051 

. 1 05108395 

. 1 05108729 

.10510905? 

.105109775 

. 1 C5 1 09555 

.105109993 

.105110297 

.10511 0595 

.106U0B9O 

.105111182 

.105111*72 

.1051 11756 

.1051 12017 

.10511231* 

.105112*69 

.105112«59 

.105113127 

.10511 3300 

.10511 3552 

.105113906 

•10511*161 

.10511**13 

. 1 0*11*656 


p 

0 

1 3 

27.-27 

18 

0,  -9 

10 

38.-36 

16 

2o,-?o 

1? 

31  .-31 

1 1 

35.-38 

1* 

?*.-?* 

10 

39,-39 

1 7 

1 3 • - 1 3 

13 

26,-?b 

12 

32,-3? 

16 

17,-17 

1 1 

36,-36 

10 

*o,-*o 

lb 

21 .-21 

1* 

28,-26 

1 3 

29,-29 

1? 

33.-33 

1 1 

37.-37 

10 

*1 .-*1 

1 1 

.38,-36 

1? 

3*. -3* 

13 

30.-30 

1* 

26,-26 

15 

22,-2? 

lb 

16,-16 

1 7 

1 * . -1  * 

18 

10,-10 

1* 

5,  -6 

20 

2.  -2 

1 1 

39,-39 

12 

36,-36 

1 3 

31.-31 

1* 

27,-27 

lb 

23.-23 

1 1 

*o , -*0 

16 

19,-10 

12 

16,-36 

1 3 

32,-3? 

1 7 

16,-15 

1* 

28.-26 

1 1 

*1 .-*1 

1? 

37, >37 

16 

2*, -2* 

1 6 

11,-11 

1 3 

37.-3* 

16 

2n,-20 

1* 

29,-29 

12 

36,-36 

1 3 

3*, -3* 

59 
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1 2 HP 

.000000121 

1 ?63 

1 ?fc>4 

.OOOUOOI 1 9 

1 ?HS 

1 2H*, 

. ooooooi 1 a 

1 ?m  r 

1 PH9 

.0000001 1 5 

12>'9 

1 ?9  0 

. 000000 1 1 1 

129] 

1 292 

.00000011 3 

129  3 

1 294 

.000000113 

1 29? 

1294 

.000000110 

1297 

1 295 

.000000109 

1299 

1 900 

.000000105 

1301 

l 302 

.000000105 

1 303 

1 3«4 

.0000001 05 

1 30  S 

1 305 

. O00000 103 

1 307 

1 30  « 

.000(1001  02 

13t>9 

1310 

. O000001  00 

lit  1 

111? 

.0  00000 '>99 

1313 

1314 

. 000O00O98 

1 3 1 S 

1 915 

.O0OU0OOU7 

131  7 

1 .3  1 F» 

.000000091 

1319 

1 32  0 

. 0 00  0UOO90 

13?1 

1 92? 

.000000090 

1323 

1 924 

.000000089 

1 32S 

1 324 

.00000008  7 

1 327 

1 925 

.000000085 

1329 

1 330 

.O00000'»85 

1331 

l 332 

. OQ00U0r34 

133  3 

1334 

.000000082 

1 33S 

l 335 

.000000080 

1337 

1335 

. 00 0 QUO  0 79 

1339 

1 3*0 

.000000079 

1 341 

13*2 

.000000077 

134  3 

l 3*4 

.000000075 

1 34-3 

13*5 

. 000000075 

1347 

1 348 

.000000074 

1349 

1 3b  0 

.000000^72 

13M 

1 3b? 

.000000071 

1 35  3 

1 3b4 

. O000O0O70 

1355 

1 3b5 

.000000059 

135  7 

I9b8 

.0000000 57 

1359 

1 360 

.000000055 

1351 

1 36? 

.000000054 

1353 

1 364 

. O0OOOO"53 

1355 

1 368 

. 000000O5? 

1357 

1 368 

.00000005] 

1 349 

1970 

.000000058 

1371 

1 9 7? 

. 0000001)58 

137  3 

1374 

.000000058 

1 375 

1374 

.000000058 

1377 

1 3 73 

,000000058 

1 379 

1360 

.000000058 

riiv,  sum,  p ^ 


. 1 08 

1 1 409  7 

1 5 

?4 

.-?5 

. 1 05 

1151 35 

1 7 

16 

.-16 

. 1 05 

1 1 5970 

1 2 

39 

.-.39 

. 1 05 

1 1 580? 

1 4 

3^ 

.-.3  0 

. 1 05 

1 1 58?« 

1 9 

7 

. -7 

. 1 05 

118054 

l 6 

21 

.-?) 

.105 

1 1 6279 

1 9 

35 

.-35 

. 1 03 

1 1 6500 

1 2 

40 

,-40 

. 1 05 

116718 

15 

26 

,-2  6 

. 1 05 

1 1 6930 

1 * 

31 

.-31 

. 1 05 

117140 

1 M 

1 2 

.-1? 

.105 

1 1 7350 

1 3 

36 

,-36 

. 1 05 

1 1 7556 

l ? 

41 

.-41 

. 1 05 

117759 

1 7 

17 

.-17 

. 1 05 

1 1 7959 

1 6 

?? 

,-P2 

.105 

1 1 81  56 

15 

27 

.-27 

. 1 05 

1 1 8252 

1 4 

3? 

.-3? 

.105 

1 1 8547 

1 3 

37 

,-37 

.105 

1 1 8729 

1 3 

38 

.-38 

. 1 05 

118909 

1 * 

33 

.-33 

. 1 05 

1 19089 

15 

28 

.-28 

. 1 05 

1 19266 

1 6 

23 

.-23 

.105 

1 1944) 

1 7 

18 

,-18 

. 1 05 

11961  1 

1 « 

13 

.-13 

.105 

119781 

1 3 

39 

.-39 

, 1 05 

1 1 9948 

1 * 

34 

.-34 

.105 

120111 

15 

29 

.-29 

. 1 05 

12027] 

19 

8 

, -8 

.105 

120430 

1 3 

40 

.-40 

' 1 05 

120589 

16 

24 

.-24 

.105 

120743 

1 4 

35 

,-35 

.105 

120895 

1 7 

19 

.-19 

.105 

121045 

15 

30 

, -30 

.105 

121193 

1 3 

*1 

.-*1 

. 1 05 

121337 

1 4 

36 

,-36 

,105 

121479 

1 6 

25 

.-25 

. 1 05 

121619 

1 8 

14 

,-14 

.105 

121756 

15 

31 

.-31 

. 1 05 

121890 

1 * 

37 

,-37 

. 1 05 

122022 

1 7 

2o 

,-20 

.105 

122151 

1 6 

26 

.-26 

. 1 05 

1 2??  77 

15 

3? 

.-32 

. 1 05 

1 22401 

1 * 

38 

. 1 05 

122523 

20 

3 

, -? 

. 1 05 

1 22640 

19 

9 

, -9 

. 1 05 

122756 

1 « 

15 

,-15 

. 1 05 

12287? 

1 7 

21 

.-21 

. 1 05 

1 22989 

1 6 

27 

,-27 

. 1 05 

1231  05 

15 

33 

.-33 

.105 

123221 

1 * 

39 

.-39 

60 

M 


? I f^f  9 V A i,  If 


ru*.  si* 


o 


9 


1 3*  I « l 

1 J«  1.  I 1*4 
1 3 * * • | Id*, 
1 3*  ),  | i** 
1 3l-  >.  1 *9i> 
1 3ti  I . 1 iv? 
) 3d  1.  1 1W4 
1 34-1  . 1 996 
1 3d  7 • 1 99  < 
1 3d  3.  1 400 

I 4o  1 * 1 407 
1 4 (J  3 • I 4 ij  4 
1*06.  UO*' 
140  7, 1*0* 

l*0-».  14  in 

I * 1 I » I 4 1 7 
1*1  3.U14 
1 * 1 6 » 1 4 1 * 
1*1  7 . I & J * 

1*16,1420 

1*21  • Ud> 
1*73.  U2. 
1*26. U26 
1*?7.14?« 

1 *?'<*•  I 4. 30 
l*ll<|4i7 
1*33. U 34 
1 4 36, 1 4 J* 
1*37.143* 
1*39.14*3 

1 *4  1 • l 447 
1 *4  3 • 14*4 
1*4*.  1444 
1447.1 4*M 
1*43.1  4b'j 

1 *S  1 . Ui)J 

1*63,  1464 
1*6*.  |4S6 
146,7. 14*5* 
1*69,1  4hn 
l*6l . 144? 

1 *f  3 . 1 4 7*4 
1*66.  Uft* 
1*67.1  46* 
1*69.14  7o 
1*71 . 1*77 
1*7  3.  1 4 74 
1*76.1  4 76 
1*77.U7H 
1*7  3. ] 4*  > 


.0vj0niK'rS4 

.oononoiS* 

. « 0 0 rt  o n i 6 1 
• noonno"* i 
. o •} n n no  6 1 
. OOOOUO<'*0 
. 000000  49 

• n . i n ,o  o o '4  x 
. ooom  O'  4f> 
.000  0 0 o I4x 
,0000  00  -44 
. OOOilOf!  f'44 

.ooooun-14  i 
. 0 o o 0 o « ' 4 ? 

. onoooo>4i 

.ooo  o 0 n i40 

.000000040 
. no o no On  3 7 
. o o n o o o o i / 

. n 0 0 o 0 f>  O 1 / 

• 0000000 3h 

.nonnooo is 

.no oo oo  34 
.OOnpono-w. 

. oonod 0 1 33 
. OOOO'IO  r 1 3 
.no ooo 003] 
.nonoonoii 
.00001)0  1 30 

• 0n000Oi>?9 
. Of. 000007* 
.000000077 
.noouono?; 

. noooon ' 7 7 
.nono  on  i'?  7 
. noooon 075 
.O0000O'l7S 
. OOOOOOi,?4 
•Oo 0000073 
.000000077 

.oonooo’.?? 

.000000 r?? 

. OUOOOO!'?! 
.OOOOOOn?! 

• noooon  ^ 7o 
.nooooOnl 9 
.000  OOO  n 19 

• o o o i)  o o n m 


. 1 OS] 7 J379 
. 1 061  734  37 
. 1 061 7 3647 
. 1 061 3 3646 
.106173747 
. 1 °6 1 7 3*46 
. 1 061  7 3944 
.1061 74040 
. 1 > 6 1 741 3? 
. 106174774 
. 1 061 7431 7 
. 1 06134400 
.1061 74466 
. 1 06174669 
.1061 34660 
. 1 06174731 
. I 061 740]  1 
. l 06174006 
.1061  74969 

. 10617603? 

. 1 061 3S1 04 
. 1 061761  73 
. 1 0617674? 
.10617*310 
. 1 061 76376 
.106176441 
.1 06136604 
. 1 061 76666 
.106176676 
.1  061 7*603 
.106176719 
. 1 06176794 
. 1 0617604* 

. 1 06)7690? 
.106176966 
. 1 0*176006 
. 1 061 76066 
. I 06)761 03 
.106176149 
.106176194 
.106176736 
.10617676? 
.106176374 
.106176366 
.106176406 
.106176444 
. 1 06)7646? 

. I 061 766)9 
.106176666 
. 1 06) 76690 


1* 

40,- 40 

16 

34,-34 

16 

26,-2* 

1 7 

2?, -72 

1* 

*1 .-41 

16 

36,-36 

1« 

16,-16 

16 

?6,-79 

16 

36,-36 

1 7 

23, .73 

1 9 

1 o , — 1 o 

16 

30,-30 

16 

37.-37 

1* 

17,-17 

1 f 

24,-24 

16 

31  .-31 

16 

36.-3S 

lb 

id, -39 

16 

3?,-3? 

1 7 

26,-75 

1H 

X 

• 

« 

X 

lb 

40,-40 

16 

11. -11 

16 

33,-33 

1 7 

?6,-?6 

lb 

41 ,-41 

16 

34,- .3* 

1» 

19,-19 

i y 

27.-77 

16 

36,-36 

70 

4.  -4 

16 

17.-12 

1« 

?n, -?o 

16 

36,-36 

1 7 

?*.-2H 

16 

37.-37 

1 r 

76,-79 

1* 

21 .-21 

16 

JM.-36 

i y 

30,-30 

16 

13.-13 

16 

J9.-39 

1H 

27,-22 

17 

31 ,-31 

16 

*0,-40 

16 

*1 .-*1 

i y 

37.-32 

16 

23,-23 

16 

1 * . — 1 * 

17 

33.-33 

.0000000)6 

.O00U00" 1 7 
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M PIRFNJ  V 6 ' 1 p CUM  # S'lM,  P J 


1 44) 

14M? 

14  4 3 

1 4 H4 

140, 

1 4*4 

1 4 4 7 

1 4MM 

1 4h  -> 

1 40'J 

145  1 

140? 

145  3 

14  94 

1 40  -> 

1494 

1 40  7 

1 40  H 

1 40  0 

140  0 

1 36  1 

1 40? 

1 SO  3 

1 404 

1 SO, 

1 404 

1 367 

1 40* 

1 30  9 

14]  0 

161  1 

141? 

1613 

1 4l4 

ISIS 

1414 

1317 

I4ln 

161  3 

1 52  0 

1 4?1 

1 52? 

1 Si?  3 

1 524 

1 3?  5 

1524 

1 6?  7 

1 420 

1 6?  5 

1 4.30 

1631 

14-3? 

1633 

1434 

163  6 

1 4 J 4 

1637 

1 4-H 

1639 

144  0 

164  1 

144? 

1643 

1 644 

154  6 

1 44S 

1547 

1 444 

1649 

1660 

1561 

1 45? 

1563 

1664 

1556 

1 454 

1547 

1554 

1559 

1460 

16M 

146? 

1563 

1 464 

1 56  6 

1 464 

1 5ft  7 

1S6-* 

1569 

] 4/6 

1 6 / 1 

15  7? 

1673 

14/4 

16/, 

14/4 

1577 

1S7S 

1579 

1 4Hn 

.onnuoo  1 1 / 
.nonooo tih 
. n n o o (>  n n i s 
•nnooon  'lb 
,nooi)0{»‘>  l s 
. no  no oo i is 
.000000014 
.OOOoOO  1 1 1 

. o n n o o o i ) j 
. nunoooni  3 
. n ij  o o o n « ] ? 
.no  no  on  n ? 
.no  non null 
. n n n o o n n 1 1 
. n o o o n o n n 
. nunooor  l o 
, n o n o o n n i o 
. 0 u 0 o 0 o n l o 
.nonooonoo 
.non noon no 
. nonoonnnq 
. nooitoonon 
. n o n n (t  n n n h 
. n n 0 o 0 0 (>  n h 
, non  non  no  7 
.OOOoOO  1 07 
.000000007 
. nuoooonnh 

.nooooonne 
, nonnnn  iOb 
.OOOOOO  n oft 
.000000005 

,000000005 

.000000005 
. 0001)00  0 05 
.000000005 
. Q0000OOO4 
,000000004 
.000000004 
. 000000004 
.000000004 
.OOOoOO 003 
.000000003 
.000000003 
.000(100003 
. 000000003 
. OOOOOO  on* 

. n o n o 0 0 n o ? 
. o o o 0 o o o n ? 
.ooooonno? 


. 1 n51 ?64?3 
. 1 051 ?6654 
. 1 041 ?44R4 
. 1 C51  ?471S 
. 1 051P5745 
. 1 051 Pb77S 

• 1 05 ] ?4B03 

. 1 051 

. 1 051 ?6R44 
.1051 26BH0 
. 1 051?6004 
. ] 051 2692M 
. 1 051?6941 
. 1 051 24973 
. 1 051 

. 1051 P7015 
. 1 05127035 
.1051 ? 7055 
. 1 051?  7n'3 
. 1 051?709? 
. 1051P71 10 
.105)-571  27 
. 1 051P7141 
,1061?715P 
. 1051?71  7? 
. 1 051  2 71  86 
,105127199 
,1  05127212 
, 1 051 ?7??4 
.105127234 
,1051?7?47 
. 1 051?  7258 
, 1 051 ?7?4R 
. 1 051 ?7?75 
, 1 051?7?«H 
. 1 0S1?7?97 
. 1 0512/305 
. 1 0 5 1 ? 73  1 4 
. 1 05] ?73?? 
. 1 051?7730 
. 1 051 ?733H 
. 1051?7345 
. 1 051? 7 351 
, 1 051 ?735  7 

• 1 05 1 ?736  3 
. 1 05127348 
. 10S1?  7373 
, 1051?/37 7 
. 1 051?  73H? 
. 1 051?738h 


1? 

? 4 

.-24 

1 / 

34 

.-34 

?0 

4 

♦ -4 

19 

14 

.“15 

1 H 

24 

.-26 

1 T 

34 

,-34 

1 7 

34 

.-36 

16 

?4 

, -26 

1 / 

37 

.-3  7 

1 o 

16 

,-16 

1 8 

2 7 

,-27 

1 / 

3w 

.-38 

1 1 

39 

.-39 

]H 

?n 

« -28 

19 

1 7 

,-l  7 

1 / 

40 

♦ -40 

IP 

?0 

,-29 

1 7 

41 

.-4  1 

?0 

4 

, - 6 

19 

IB 

,-]  8 

1 8 

30 

« -30 

1« 

31 

« -3 1 

19 

19 

.-19 

18 

3? 

.-32 

1 8 

33 

,-13 

1* 

?0 

♦ -20 

IB 

34 

.-34 

?0 

7 

, -7 

19 

2l 

.-21 

1 8 

36 

,-36 

18 

34 

, -3b 

19 

?? 

.-2? 

18 

37 

,-37 

16 

38 

.-38 

10 

?3 

,-23 

1 8 

39 

,-39 

? o 

M 

. -8 

1 B 

4 0 

,-40 

19 

2 4 

,-24 

1R 

41 

,-41 

19 

25 

,-25 

10 

26 

,-26 

?n 

Q 

, -9 

19 

? 7 

.-27 

19 

?R 

,-2R 

10 

29 

,-29 

?() 

In 

, -1  0 

1 0 

30 

.-30 

1 o 

31 

,-31 

19 

3? 

,-32 

62 


0U*.  S'.IM 


N 

1 “»M  • 1 *M2 
I SP  ! • 1 SH« 
1 5**  3.  ) SH* 

I S*7.  I **M 
1 Sh-9.  m^u 

1591  . 1*92 
1 So  4, 1 s^4 

ISO-,,  I *9s 
ISO  7 , 1 *'9M 
1 599.  1 *,00 
ISO  I .1*0? 
ISO  3.  1 *,04 
1 60S,  1 *0* 
1*07. 1 *,  0 h 
lbOOt  1*1(1 
1 «>  1 1 .1*1? 

1 ft  1 1, 1*14 
IMS.  1*1* 
1M  7.  |*1* 

1 h ) 0.1*20 

lb?l . 1*2? 
lb?  4,  j*2* 
lb?S. 1*2* 
1*2  7. 1*?9 
1*29.  |*Jf. 
1*31 . 1*32 
1*33. 1*34 
ISIS. 1*1* 

1 b 1 7 , |*1* 

1 b 1 9 , l**i> 

1 b<*  1 .1*42 
1*41, 1*44 
1 *4S, 1 *4* 
1047,1*4* 

1 S4  9. 1 *SO 
1*M  » 1*22 
lb*3. 1*24 

l b*3. 1 *3* 

1 bs  7 , 1 *s* 

1 **9, 1 *60 


?I^FO  W/k.'tF 

.0  0 o<)  0000? 

. OOOuUO'tO? 
. O000000O? 
.000000002 
. o o o o n o o o i 
.ooouonooi 
. o 0 0 0 0 0 0 o 1 
. 0 o 0 0 o 0 0 0 1 
.000 n oo 001 

. OOO000O01 

• 0 0 0 1)  0 0 o 0 1 
. 0000 0 0 0 0 1 
. 0 0 n o 0 0 O 0 1 
.000000')  01 

• o U 0 0 0 0 0 0 1 
.OOooyooOl 
.OOOOOOoOO 
. OftOoOOdOO 
.OOOOOOOOO 
.OOOOOOOOO 
.OOOOOOOOO 

.oooooonno 
. ooooooooo 

.O()0u00000 
.OOOOOOOOO 
.Ooooooooo 
. ooooooooo 

.OOOOoOoOO 

.OOOOOOOOO 

.OOOOOOOOO 

•OOOOOOOOO 

.OOOOOOOOO 

. ooooooooo 

.OOOOOO'lOO 

.ooooooooo 

.ooooooi-no 

.OOOOOOOOO 
.OOOOOOOOO 
•OOOuOO^OO 
. OOOOOOOOO 


. 1 OS!  27  389 

• 1 OS! ? 7191 
.10*12739* 
. I 0S1273Q9 

• 1 0*1 274  0? 
.10S12740S 
. 1 0*1 274  0P 

• 1 0S1  ?74l(l 
. 1 0*1 ?74l 3 
.10Sl?74lS 
. 1 0S1 2741 7 

• 10*12741*9 
. 1 usi 274?1 
. 1 0*1 274? 1 

• 1 nS1274?4 
. 1 0S1  274?S 
.10*1274?* 
,10*1274?7 
. 1 0*1 2 7 4 ? P 
.10*1274?9 
. 1 0*1  2 7429 
. 1 0*1  27430 
. 1 0*127430 
. 1 0*1  27430 
.10*12743] 

. 1 051 27431 
. ! 0*127431 
. 1 0*127432 
.10*12743? 
.10*1?743? 
.10*12743? 

. 10*127433 
. 1 0*127433 
. 1 05127433 
.105127413 
.10*127433 
. 1 0*1 27433 
.10*127434 

. 1 05127434 
. 1 0*1 27414 
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?0 

11. -11 

19 

33.-33 

19 

34 . -34 

19 

35.-35 

20 

12.-12 

19 

3*. -3b 

19 

.37,-37 

19 

3® . -38 

?0 

13,-13 

19 

39,-39 

19 

40.-40 

19 

*1 .-41 

20 

14,-14 

?n 

is. -1* 

?o 

l*,-lb 

20 

17.-17 

?0 

19.-18 

20 

13.-19 

?0 

?0,-?0 

20 

21 .-21 

20 

2?. -22 

20 

23,-23 

20 

24,-?4 

?0 

2*.-?* 

20 

2*. -2b 

?0 

27, -?7 

20 

29,-28 

20 

29,-29 

20 

30.-30 

20 

31,-31 

20 

3?. -.3? 

20 

33,-33 

?0 

34,-34 

20 

3*. -35 

20 

3*. -3b 

20 

37,-37 

20 

3". -38 

20 

39,-39 

20 

SO.-40 

20 

*1 .-41 
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Table  IV 

CENSORING  System  Probabilities 


tlie  probability  that  an  aperture  of  diameter  D/r0 
will,  during  a single  short  exposure,  receive  a wavefront 
whose  rm  s distortion  over  the  aperture  (with  tilt  not  con- 
sidered a form  of  distortion)  will  be  less  than  one  radian. 


»/r0 

P 

1 CENSOR 

2 

0.  986  ± 0.  006 

3 

0.  765  ± 0.  005 

4 

0.  334  ± 0.  014 

5 

(9.  38  ± 0.  33)  x 

io- 2 

6 

(1.  915  ± 0.  084) 

x lO-2 

7 

(2.  87  ± 0.  57)  x 

10-3 

10 

(1.  07  ± 0.  48)  x 

10-8 

IS 

(3.  40  ± 0.  59)  x 

io-  ‘5 
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5 Rome  Air  Development  Center  ? 


RADC  is  the  principal  AFSC  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences , intelli- 
gence, command,  control  and  communications  technology, 
products  and  services  oriented  to  the  needs  of  the  USA.F. 
Primary  RADC  mission  areas  are  communications , electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects,  intelligence  data  collection  and 
handling,  information  system  technology,  and  electronic 
reliability,  maintainability  and  compatibility.  RADC 
has  mission  responsibility  as  assigned  by  AFSC  tor  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence , mapping,  charting,  command, 
control  and  communications  areas. 
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